{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# GNN Explainer\n",
    "\n",
    "This notebook is designed to visualize the results of the GNN Explainer."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import os\n",
    "import networkx as nx\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Configuring the experiment you want to visualize. These values should match the configuration:\n",
    "\n",
    "> TODO: Unify configuration of experiments in yaml"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [],
   "source": [
    "logdir = '../log/'\n",
    "expdir = 'syn2_base_h20_o20_explain'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the produced masks"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['events.out.tfevents.1571881971.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571933062.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571878932.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571932752.Dylans-MacBook-Pro-2.local',\n",
       " 'masked_adjsyn2_base_h20_o20_explain.npy',\n",
       " 'events.out.tfevents.1571878946.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571879199.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571879111.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571881551.Dylans-MacBook-Pro-2.local',\n",
       " 'masked_adj.npy',\n",
       " 'events.out.tfevents.1571878892.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571878623.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571881440.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571882109.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571882039.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571881076.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571932827.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571878863.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571878915.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571881416.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571878403.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571932956.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571876351.Dylans-MacBook-Pro-2.local',\n",
       " 'events.out.tfevents.1571932775.Dylans-MacBook-Pro-2.local']"
      ]
     },
     "execution_count": 53,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dirs = os.listdir(os.path.join(logdir, expdir))\n",
    "dirs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "masked_adjsyn2_base_h20_o20_explain.npy\n",
      "masked_adj.npy\n"
     ]
    }
   ],
   "source": [
    "masks = []\n",
    "# This would print all the files and directories\n",
    "for file in dirs:\n",
    "    if file.split('.')[-1] == 'npy':\n",
    "        print(file)\n",
    "        masks.append(file)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [],
   "source": [
    "def show_adjacency_full(mask, ax=None):\n",
    "    adj = np.load(os.path.join(logdir, expdir, mask), allow_pickle=True)\n",
    "    if ax is None:\n",
    "        plt.figure()\n",
    "        plt.imshow(adj);\n",
    "    else:\n",
    "        ax.imshow(adj)\n",
    "    return adj"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 76,
   "metadata": {},
   "outputs": [],
   "source": [
    "def filter_adj(adj):\n",
    "    filt_adj = adj.copy()\n",
    "    filt_adj[adj<0.8] = 0\n",
    "    return filt_adj"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 91,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA20AAAEyCAYAAABptTjBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdeVhU59k/8O/AIOCCihvIIiggLiBbImnzJqY2aZukagwRd9xC9sa8TRrTvG9r2xhj2jRpWlN/vHEBFZSYqFmMJpqaJk1ccN8ZEBQQjaioKNsM5/fHMIcZmH09c+b7ua5chZlznvPMYJ9r7rnvcz8KQRAEEBERERERkST5eXoCREREREREZBqDNiIiIiIiIglj0EZERERERCRhDNqIiIiIiIgkjEEbERERERGRhDFoIyIiIiIikjAGbUREREREbjBu3Di8//77np4GAGDNmjW4++67PT0NshKDNiIiIiIiG3z77bf40Y9+hN69eyM0NBQ//vGPsX//fk9Pi2RM6ekJEBERERF5ixs3buDhhx/GP//5T0yZMgUtLS345ptvEBgY6NZ5qNVqKJX8KO8rmGkjIiIiIrJSaWkpAGDatGnw9/dHcHAwHnjgASQnJ2Px4sWYOXOmeGxlZSUUCgXUarX4WHl5Oe68806EhIRg4sSJuHr1qvhcQUEBhgwZgn79+uFPf/oTYmJisHPnTgDA4sWLkZWVhZkzZyIkJARr1qzBvn37cNddd6FPnz4IDw/Hs88+i5aWFnE8hUKBd999F0OHDkX//v3x0ksvoa2tzeD1vPjii+jbty9iY2Px+eefu+Q9I8cxaCMiIiIislJCQgL8/f2Rk5ODzz//HNeuXbPp/IKCAqxatQq1tbVQKpX41a9+BQA4efIknn76aaxfvx61tbW4fv06ampqDM7dunUrsrKyUF9fjxkzZsDf3x9vv/026urq8P3332PXrl147733DM7ZvHkzSkpKcPDgQWzduhWrVq0Sn9u7dy+GDx+Ouro6/OY3v8H8+fMhCIKd7wy5EoM2IiIiIiIrhYSE4Ntvv4VCocDjjz+OAQMGYMKECbh06ZJV58+aNQujR49Gjx498Kc//QnFxcXQaDTYtGkTfvnLX+Luu+9Gt27d8Mc//hEKhcLg3LvuuguTJk2Cn58fgoODkZ6ejszMTCiVSsTExOCJJ57A119/bXDOyy+/jNDQUERHR2PhwoUoKioSnxsyZAgef/xxMQitra21+nWQezFoIyIiIiKywYgRI7BmzRpUV1fj+PHjuHDhAhYuXGjVuVFRUeLPQ4YMQWtrK+rq6nDhwgWD57p3745+/fqZPBfQlmo+/PDDCAsLQ0hICH7729+irq7O7PUuXLgg/h4WFmZwPQBoaGiw6nWQezFoIyIiIiKyU2JiIubMmYPjx4+jR48euH37tvjcxYsXuxxfVVUl/nz+/HkEBASgf//+CA8PR3V1tfhcY2Mjrly5YnBu58zbU089hcTERKhUKty4cQOvv/56l/LGztcbPHiwfS+UPIpBGxERERGRlU6fPo233npLDLCqqqpQVFSEzMxMpKSk4N///jfOnz+P69evY+nSpV3OX7duHU6ePInbt2/jd7/7HbKysuDv74+srCx88skn+O6779DS0oLFixdbvL/s5s2bCAkJQc+ePXH69Gn885//7HLMn//8Z1y7dg1VVVX429/+huzsbOe8EeRWDNqIiIiIiKzUq1cv7N27F2PHjkWPHj2QmZmJ0aNH46233sL999+P7OxsJCcnIz09HQ8//HCX82fNmoU5c+YgLCwMTU1NePfddwEAo0aNwt///ndMnToV4eHh6NmzJwYOHGh2K4G//OUvKCwsRK9evfD4448bDcgmTpyI9PR0pKSk4KGHHsL8+fOd92aQ2ygEtoghIiIiIpKUhoYG9OnTByqVCrGxsXaNoVAooFKpEBcX5+TZkbsx00ZEREREJAGffPIJbt++jVu3buHFF19EUlISYmJiPD0tkgAGbUREREREErB161YMHjwYgwcPhkqlwoYNG7o0HyHfxKDNSy1evBgzZ84EAFRWVkKhUECtVjtl7DVr1uDuu+8Wf+/ZsyfOnj3rlLHlTv/vQuSLuDZJE9cmIu/w/vvvo76+HtevX8euXbswfPhwh8YTBIGlkTLBoM3DYmJiEBwcjJ49e4r/6e+f4Uzjxo1D37590dzcbNN5DQ0NGDp0qEvmRETSxLWJiIhIOhi0ScAnn3yChoYG8T9X7J9RWVmJb775BgqFAh9//LHTxyci+eHaREREJA0M2iRq9+7diIyMNHgsJiYGO3futGu8goICZGZmYs6cOcjPzzd47sqVK5gwYQJCQkJw5513ory83OB5hUKBsrIyAMBnn32G1NRUhISEICoqCosXLzY49ttvv8WPfvQj9OnTB1FRUVizZg0AoLm5GS+++CKio6MxaNAgPPnkk2hsbDR4rW+99RYGDhyI8PBwrF69WhyzsbERv/71rzFkyBD07t0bd999NxobG/HQQw/h73//u8H1k5OTsXnz5i6vX1emtXr1akRFRaFv375YsWIF9u/fj+TkZPTp0wfPPvuseHx5eTl+8pOfoF+/fujfvz9mzJiB+vp68flly5YhIiICvXr1wvDhw7Fr164u12xtbcW0adPw6KOPoqWlxdSfhsircG3i2kRERO7HoM1HFBQUYMaMGZgxYwZ27NiBS5cuic8988wzCAoKQm1tLVatWoVVq1aZHKdHjx4oKChAfX09PvvsM/zzn//Eli1bAADnzp3DL37xCzz33HO4fPkyDh8+jJSUFADAokWLUFpaisOHD6OsrAw1NTX44x//KI578eJFXL9+HTU1NVi5ciWeeeYZXLt2DQDw4osv4sCBA/juu+9w9epVvPnmm/Dz80NOTg7WrVsnjnHkyBHU1NTgoYceMjn/vXv3QqVSYePGjVi4cCGWLFmCnTt34sSJEyguLsbXX38NQFsD/sorr+DChQs4deoUqqqqxA+BZ86cwT/+8Q/s378fN2/exI4dO7p0dmpsbMSkSZMQGBiI4uJidOvWzYq/EpHv4dqkxbWJiIjMEsijhgwZIvTo0UPo3bu30Lt3b2HixImCIAjCv/71LyEiIqLLsV9++aUgCILw+9//XpgxY4YgCIJQUVEhABBaW1uNXuObb74RlEqlcPnyZUEQBGH48OHCX//6V0EQBEGtVgtKpVI4deqUePwrr7wi/PjHPxZ/ByCoVCqjYz///PPCwoULBUEQhNdff12YNGlSl2Pa2tqE7t27C2VlZeJj3333nRATEyO+1qCgIIP5DxgwQPj+++8FjUYjBAUFCYcPH+4ybmNjo9CnTx+htLRUEARB+PWvfy089dRTRuepe4+qq6vFx0JDQ4UNGzaIv0+ePFl4++23jZ6/efNmISUlRRAEQVCpVMKAAQOEL7/8UmhpaTE47ve//73wy1/+UrjnnnuE5557TmhrazM6HpHUcW3i2kRERNLBTJsEbNmyBfX19aivrxe/GXam/Px8PPDAA+jfvz8AYPr06WIZ0uXLl6FWqxEVFSUeP2TIEJNj7d27F/fddx8GDBiA3r17Y8WKFairqwMAVFVVYdiwYV3OuXz5Mm7fvo309HT06dMHffr0wc9//nNcvnxZPKZfv35QKpXi7927d0dDQwPq6urQ1NRkdNygoCBkZ2dj3bp1aGtrQ1FREWbNmmX2vRg0aJD4c3BwcJffGxoaAACXLl3C1KlTERERgZCQEMycOVN8nXFxcXjnnXewePFiDBw4EFOnTjVo0LBnzx4cPXoUixYtYpte8mpcm7g2ERGRNDBok6gePXrg9u3b4u8ajcbgg4S1GhsbxdKasLAwhIWF4e2338aRI0dw5MgRDBgwAEqlElVVVeI558+fNzne9OnTMWHCBFRVVeH69et48sknIQgCACAqKqrLPScA0L9/fwQHB+PEiRPiB8Dr16+LH0LM6d+/P4KCgoyOCwA5OTlYv349du3ahe7du+Ouu+6yOKY1fvvb30KhUODYsWO4ceMG1q1bJ75OQPs+fPvttzh37hwUCgVefvll8bkHHngAr7zyCsaPH29Q6kUkB1ybOs7l2kRERO7CoE2iEhIS0NTUhM8++wytra147bXXbG6HDWi/Kff398fJkydx+PBhHD58GKdOncJ//dd/oaCgAP7+/pg8eTIWL16M27dv4+TJk12aAei7efMmQkNDERQUhH379qGwsFB8bsaMGdi5cyeKi4uhVqtx5coVHD58GH5+fnj88cfxwgsv4IcffgAA1NTUYMeOHRbn7+fnh3nz5uG///u/ceHCBWg0Gnz//ffie3HXXXfBz88Pv/71ry1+k22LmzdvomfPnujduzdqamrw5z//WXzuzJkz+Oqrr9Dc3IygoCAEBwfDz8/w/0q/+c1vMH36dIwfP178FpxIDrg2aXFtIiIid2LQJlG9e/fGe++9hwULFiAiIgI9evTo0rHNGvn5+Zg7dy6io6PFb7PDwsLw7LPPYv369VCr1fjHP/6BhoYGhIWFYc6cOZg7d67J8d577z387ne/Q69evfDHP/4RU6ZMEZ+Ljo7Gtm3b8NZbbyE0NBQpKSk4cuQIAG1Hs7i4OGRmZiIkJAQ//elPcebMGatew1/+8hckJSXhjjvuQGhoKF5++WW0tbWJz8+ePRvHjh1z6saxv//973Hw4EH07t0bDz30ECZPniw+19zcjEWLFqF///4ICwvDDz/8gKVLl3YZ43//938xadIk/PSnP8XVq1edNjciT+La1IFrExERuYtC0K+rIOqkra0N/v7+OHfuHKKjoz09HaMKCgqQl5eHb7/91tNTISI34dpERES+RGn5EPJlx48fR1BQEMLCwjw9FaNu376N9957D08//bSnp0JEbsS1icg96hqaselANU5fvIEbTWqEBCmRGBaCx9Ij0a9noKenR+QzGLSRSR9++CFyc3OxbNkySe7ls2PHDkyePBk//elPMX36dE9Ph4jchGsTkesdqarH8t1l+LpU22ioWd1R+hukvIi3d5Zi3PABePreOIyJ6uOpaRL5DJeVR27fvh3PP/88NBoNFixYgEWLFrniMkRENuHaRERk3ro9lViy7TSa1BqY+5SoUABBSn+8+mAiZmbGuG1+RL7IJUGbRqNBQkICvvzyS0RGRuKOO+5AUVERRo4c6exLERFZjWsTEZF52oDtFBpb2ywf3C44wA+vPjiCgRuRC7mke+S+ffsQFxeHoUOHolu3bpg6dSq2bt3qiksREVmNaxMRkWlHquqxZNtpmwI2AGhsbcOSbadxtLreRTMjIpfc01ZTU4OoqCjx98jISOzdu9fgmLy8POTl5QEADh84jO4IMXi+ZWhwR0q+TYHAc9rNXJtjg8VjAisau5wDAN3OGj5ui9ZhQdpLtinE8Ztjg7tcy9s0x3QHAARW3hbfQ2teU8vQYIfeT3dqjumOwMrbFo9zJWf8G5SLJtxCi2D7/l2uZM3aBFhen4jIu0lxfZKC5bvL0KTW2HVuk1qD93aXYcXMDCfPiogADzYiyc3NRW5uLgAgRBGKsYrxBs+Xv5YKf6V24WjT+GHo9MMAgJq/jkLE5BPaY4pSMKz98YqiMQgObgEAhE06ZfU8zhamiGNrH9B7UtH+v5UdP3c53luca/9fBbSvR/ezBeVLOt5jQPs+A0DstCPOnJ1znINVr8ladv2tK7T/o1qbBgCIn33QeRPyMnuFXZ6egt0srU9E5N28eX1ylbqGZnxdetnsPWzmCALwrzOXcaWhmV0liVzAJUFbREQEqqqqxN+rq6sRERFh0xiamwFQhGhXDo3aD+eKkwAAQyYfw4XN2vtP/Jo6VpbYaUdQti4VAFC5MRkx2UcBAOWFKQCAXj2acL1em22KzzmIi1tGAACGTur4UH6uOAltGj9xPGMUio5rlq3VXi+4R0tHILk+FcNmHLL4+lQFHR/qKzcmA4A457K1qYib1TGG7vUOfuRkx/n5aVD4aeeif2xnutevC7w6ByKqfO08AoLU4vV1xwGAptlfPEZQ+0EpqLXnrUlHt+7aIFmt9teOEaBB9GPHzL/u/DQE9dCe19QQiIDgVoPXrqN7T1qbtP9E43NMBz/GAsnzHySJczlbmII2jTaaM/deAdr3q02t/TcQP/0gSv/fHQAA/55q8e+qe28AGLyX+u+tv9K20hJyD2esTUREcrTpQLXDYygAbDpYjSfuGeb4hIjIgEvuabvjjjugUqlQUVGBlpYWbNiwARMmTHDFpYiIrMa1iYjIuNMXbxi09bdHk7oNp2tvOmlGRKTPZS3/t23bhoULF0Kj0WDevHl49dVXTR4bNCwCkUufRJvGH5qbAQCAhCf34fwH2uxaa6u/QdWbfnZDP0uln3HRz2Tp6DJjnbMtuvPULf5iRqdzRkp/7KpNowEAUVnHAQClq9ORMPeAeKyYFZzSkXUqfT8DCQtKAACXPx6OARPOmHw/7CHpskUApSvuBKD9u+rTZa1MlSGWr9f+zazJXpqjy6zaUjqrr3TFnV3mTubtFXbhhnDV09Powpa1CWB5JJEcSXV98qR5+fvx1ekfHB5nfOJArMy5wwkzIiJ9LgvabBE4NAIRrz8Nf/82sTxRqVdqV7kxGQEB2vvbNBo/dAvQlug1NAQhoJv2Z4UC4vEVRWMAheHLip16VAym1K1KgyBAV1YZN9N4YKBfDqe0ogSwdKX2JtyE+SUmj6nYkCzOu3NpoP7zsVO7Ptflev93BxIe32/0uc6llZ1LL43OrWgMAoNaDc4zpnNZp/5jGo2fwb1w1jA2nrt0LiO1d07GviyQE2vLfwH5fChi0EYkP3JZn5xp4cZD2HL4gsPjPJISgbezUywfSEQ2cUl5JBERERF5j8SwEAQqHftYGKT0Q2J4LyfNiIj0SSPTFhOJsN89B2V3NTTtTSD8lG3w99fWVncufbTEWDbgwuaRZrNGgDaromnP9GkaApCQq81elRem2Jw10qn+cBQAIPLRE0afrygaY/Y1qfLTDJpw2JKBk7rywhR0765tuXyrIQjdArVZU0uZTKCjgUp8zkGUrk4HAIMSVUfoj+3wWAVpss26WSKXb7KZaSOSH7msT85U19CMHy/7yqH72gKVfvju5Z+weySRC3is5b++wHO3xfu9dCWMGrW/WBKpH9jo/1y2LhWK9pvd9IM0ZYBavH9NaNMeEP+I4QdnY2VspkrgLAVstpSM6YiBxjTzgYayPZDRMRas6XfLNMVU2Z5+gGJLYNw5mLSk5iNt8KrrsgkYvq+l/3cHotvLU89/kARNe7dHU8Gp/rVNBWvG7i00pkvJqGB534DSVdoSWL8AjdlyU18N2IiIyLv07xmIexMG4MtTl+xq+69QAPcNH8CAjchFWB5JRERERHhmXByClP52nRuk9MfT4+KcPCMi0pFEpq05Nhjlr6cgKKgVQyZ37K2l0dszTT8LpGsOEZd9SNxHC+goHdQ0+0NoPzd+jvEsjC4DUv3hKLQ0aztW6ncwPP9BEgLaG56EW+g42DnLplqjzaL5BbQBTcbPMZYd0nVK7NWzEQ23grRzmqq3p9qadKB9bzb9DI6ueYvu9RgrxdQ/XteVM/qxYwYZK2PZzM6vS/d+WpNl02X3hFY/JEw23ZQFAPyC1R1NPx7rml3T7+ZpKrPZuStny+1uXY4pW5fapeGMwt/wK0Xda+x8rH4ny4R55l+PTkXRGHGvOYW/wMwbERFJ1pioPnj1wUQs2XYKja3Wl0kGB/jh1QcTkRzZx4WzI/Jt0rinbWgEIpY8g2EzDhl0OzTWpl2/FNBgE20TH+R1nSG792hGa6v22yNL5XLG6MotuwWpu9xz1fnapjoR6n/ot1SKaKycUP/1dAvsOg+z87fQIdOV7Ckf9RRznTh1KjYk23RPodS3Y3AVudwzwnvaiORHLuuTq6zbU4kl206jSa0xWyoptLUhuJsS//PQCMzMjHHb/Ih8EcsjiYiIiEg0MzMGG3Mz8bORgxCo9ENQp66SQUo/BCr90L+pBj9VHGPARuQGksi06b7JLi9MgV/7/mqx044YZIf0fza2d5alLowkbZ03KCfvJ5dvsplpI5IfuaxP7nCloRmbDlbjdO1N3GhqRUhQABLDeyErLRI3Ll9ARkYGTp06hYEDB3p6qkSyJqmgTV/lxmQo2gO4IVOOiaWFmgYlEp4wX75mSXlhCnT9AYea6QxpaWPl0vfbN9FeUGJTCVzNR6O6lD0SyY1cPhQxaCOSH7msT1Lw/PPPQxAEvPvuu56eCpGssTySiIiIiOzy6quvYv369Th79qynp0Ika5LKtHUucdTvcqhPl3Xz82uzqSTSmo2pa9ubn+h3jKwoGgN/ZZvRuZij3/HQXc5/kGR2jsa6J9rC3MbTnni9JF1y+SabmTYi+ZHL+iQVixcvRnl5OdauXevpqRDJliSCtqBhEYh8/SkEB7egoUHb6t7f33hAZmvnPksqisYAuvvo9MY9V5xktMtk+fpUtLVv/Mz27VriRuEuvCdNPyB09r8BuXPH38cYuXwoYtBGJD9yWZ+k4ubNm4iPj8cXX3yB5ORkT0+HSJZYHklEREREduvVqxd++9vf4pVXXvH0VIhkSxKZNv1vsvWbf+iae7S1+UHZvtF1a7MScbMc3/PrXLG29FI/m1a6MgMJ8zs2TbbUXMRU+aYxqvw0g5LC0lXtTUys3KQZMNyXTl95YUqXPeFMqdo0Go3XtdlM/dcqdfobe7uapzJTciOXb7KZaSOSH7msT1LS3NyMxMRE5Ofn45577vH0dIhkR3JBm46qIA3KbtpALXbqUXFza0Hj57YP7+6mH0iau3dMx5p79KhrME7uIZcPRQzaiORHLuuT1Kxbtw7Lly/Hd999B4VCYfkEIrIayyOJiIiIyGHTpk3DrVu3sHXrVk9PhUh2lJ6eQGflhSkAgPjpB8VMEgAIbe3NP+YcMNho2xOcneHSdcMcNqXj9egybOY6MtpyfWsyd55W89EoALBqDzvdvxNrykIT5pcwK0lERORi/v7+WLp0KV566SX88pe/hL+/v6enRCQbzLQRERERkVM8+OCD6N+/PwoKCjw9FSJZkVymrVePJvFnXVbkwuaRiH+kI0PUvUez2+elz95sjanskOa26T+DI/uenf8gCU03AwEACTnSv6fLUoZNvzGMtY1XdJhhIyIicj2FQoFly5YhOzsb06ZNQ1BQkKenRCQLksu0Xa/vjuv13VG2NhXnipNwrjgJgx85aXBMa6s/WlsNU+76pZSW6IInoKMDpDln9Y43RtfxsvPP+srWpWLY9MNGg42Ex/cj4fH94u+6MkEAKH0/w+L8TIl+7BgS5pXY1KHSEt3fxJTy9aliuae9dKWcnbXeCkDrrQCDxyz9bazljHkbY+69kitdcE1ERL7prrvuQlpaGpYvX+7pqRDJhuSCNiIiIiLybkuWLMGyZctw/fp1T0+FSBbsbvlfVVWF2bNn49KlS1AoFMjNzcXzzz+Pq1evIjs7G5WVlYiJiUFxcTH69u1rdixTLbXFBh0zTDccsaUhhaVjzTX9IHmr/lCb3Yx81HITFLKOJ1tqu2N9IiLvxZb/7jF37lxERETgtdde8/RUiLye3UFbbW0tamtrkZaWhps3byI9PR1btmzBmjVrEBoaikWLFuGNN97AtWvXsGzZMrNj6X8ourhlBAAgbNIpg2NUBe3dD2cbdj/UlcfpB1vnipMMNs02p6JojMnNs00xdk1HXNg8EgC6lIG6Q9Wm0QCAqKzjbj2X5M+TH4pctT4RkTwwaHOP8+fPIzU1FcePH0dAr1BsOlCN0xdv4EaTGiFBSiSGheCx9Ej06xno6akCAOoamiU/R/JddjciCQ8PR3h4OACgV69eGDFiBGpqarB161bs3r0bAJCTk4Nx48ZZ/FBERORMXJ+IiDwvOjoaE+Y+h0l/+Qz13SMAAM3qNvH5IOVFvL2zFOOGD8DT98ZhTFQfj8zzSFU9lu8uw9ellyU7RyK7M236Kisrcc899+D48eOIjo5GfX09AEAQBPTt21f83RTdN9nWlCdWfzjKZAmbfndBY3RNR6IfO2bxWG9V89Eoq/Y581WqNemIn3PA09PwCVL5JttZ6xMRyYdU1ie5W7enEq99dgqNLWoo/Ey3UVAogCClP159MBEzM2PcN0Fo57hk22k0qTUw94nYk3MkApzQ8r+hoQGPPvoo3nnnHYSEhBg8p1AooFAojJ6Xl5eHvLw8AEArtC389QO2iqIx0Gi0/wfX30S7pdmwe6DhBc3Hn9GPdZRMWhOsWdqQ2dllktbcV2Up2PRkwKYqSOtSvmqJrruiteWs9tAvlzUVsJWtbd+wfZZnNmwn13Dm+kRERNbTBkOn0KRuMxuwAYAgAI2tGizZpr01xl1BkW6Oja1tFo/11ByJdBzqHtna2opHH30UM2bMwOTJkwEAgwYNQm1tLQDtfSUDBw40em5ubi5KSkpQUlKCALBOmIici+sTEZFnHKmqx5Jtp60KhvQ1trZhybbTOFptvgLCGbxhjkT67C6PFAQBOTk5CA0NxTvvvCM+/tJLL6Ffv37ijf5Xr17Fm2++aXYs/fIjXealpTEA8TnarE3lxmTEZFu/OXLpygwkzDe9N5mvd4mUa2koSYsny49ctT4RkTywPNK1cteW4MtTl8yWG5qiUAA/GzkIK2bav0+tNbxhjkT67M60/ec//8HatWvx1VdfISUlBSkpKdi2bRsWLVqEL7/8EvHx8di5cycWLVpkcazWYUFiaWCbxg9tGj8xYAMglkmaU9vedRKA2YANcF45ozN4YvPl2GlHPBaw6f+ddB1BOztbmIKLW0aInUSdrWJDsk2bsZP3ceb6RERE1qtraMbXpZftCoYAbRniv85cxpUG15Wme8MciTqz+562u+++G6aSdLt27bJ7QkREjuL6RETkGZsOVDs8hgLApoPVeOKeYY5PyAhvmCNRZw43InGGtjYFmhq7Aego2dMvYdQ0dDQfOf9BkkFDEZ3wSafc0tTCEltLLz05V1PKC1PEDcgd6UZprLFKuN7+e6aalrg6E2qqsYyrVWxI9ti1iYiI3OH0xRsGLfPt0aRuwyf/Pp58+i4AACAASURBVAC/M1/Bz8/Pof/8/f27PPbt8VtOmePp2psOjUFkC0kEbYEVjQYdIgHDD+4JufvFnwMC1AbHVW7UlrnFZB8VAyBbN8y21CXSFu4qvTS39YE+U5uSmzNM7zU40o3SmvnZw54N0aWAARsREcndjSa15YOscOnaTew7tw9tbW12/6fRaIw+fmXUY0D/BIfneKOp1QmvlMg6kgjaiIiIiMj7hQQ556Pl3Xem4+3s+U4Zq7OFGw9hy+ELDo8TEmRmGyoiJ3Oo5b+zNMcGi3ue6TtbmIKzhSko13tOv7wO0GbYYrKPih0RAdNdEXXjdRY79ahLsiC6LKD+750fs4Uqv6Nxh7ksVumKO1G64k4A2gybrXunSZ03ZtmIiIh8QWJYCAKVjn28DFL6ITG8l5Nm1JU3zJGoM0kEbYEVjV3KCvWDsGGdnitbmypuhqzjr7Rcmzx0+mG3do6MyT6Kc8VJ4r12ugCzatNou8bT76hpTsKT+5Dw5D6jz5WvTzX6uFyc/yAJ5z+wvyOnLtglIiIi22WlRzo8hgAgK83xcUzxhjkSdSaJoI2IiIiIvF//noG4N2EAFHaer1AA9w0fgH49A506L33iHO2cpDvmSNSZpO5pO1uYAoVC26Y7dtoRkxmTbkFdb3KNfuyYmEUaNuNQl+fdpfNG4Ma6Q0ZlHe9yDgCbNhC3lyveG2c2cnGUsc6itjCVoSQiIiLrRF4/gTZ1LyiUtgc1QUp/PD0uzgWzMvTMuDh8o6pDY6vG5nPdNUcifZIK2vRLF6s2jUZ0lvEP4KY+mLdpun5lYiqg0N3b5uxySVsDr+oPRyHmUc8HO46QQrBGREREnqVWq/Hiiy/i888/x6/eXIv3D1xBY6v1rfWDA/zw6oOJSI7s48JZao2J6oNXH0zEkm2nJDtHIn2SCtqIiIiIyPvU19cjOzsbALBnzx707dsXgwZVYsm202hSayAIps9VKLTZq1cfTMTMzBj3TBgQryXlORLpKATB3D9R9whRhGKsYjwAiA1G4mZZV8ZX+n4GACBhQYn4mKkNuE1xVdbNW5WvT/VoiamUSGHDdm+1V9iFG8JVT0/DYfrrExHJg1zWJ6koLS3FhAkT8LOf/QxvvfUWlMqOnMDR6nq8t7sM/zpzGQpoN6XWCVL6QYD2/rCnx8V5LHvlDXMkYqaNiIiIiOzy5ZdfYubMmfjTn/6E3NzcLs8nR/bBipkZuNLQjE0Hq3G69iZuNLUiJCgAieG9kJUW6fGGHt4wRyLJZdpqPhoFALh1PQgJcw+Ix+g3GbG34Yg7G36Qc7giC6prcONo0xIyTy7fZDPTRiQ/clmfPEkQBCxfvhyvvfYaNm7ciHvvvdfTUyKSNckFbbbS7efmrRsuny1MEQMSSwFKxYZkNv0gryGXD0UM2ojkRy7rk6e0trbiueeew7fffotPPvkEsbGxnp4SkeyxPJKIiIiIrFJXV4esrCz06tUL3333HUJCQjw9JSKfILnNtcvXp4rlj9aInXbE6ixb5cZksUTSEl3Wy9gYjirXG1s/qxbQTY2Abl33oNMxl2W7uGUELm4Z4fDc7HG2MMXk+2Uv1Zp0qNaki7/Xmnlt54qTxIYhxlzYPNKpcyMiIvJFJ06cwNixY5GZmYktW7YwYCNyI8mWR54rTkLzTe1NnwnzOzpDqtakw6+bdiPEYRLv9qhf+kjeQbUmHfFzDlg+kCySS/kRyyOJ5Ecu65M7ffrpp5g3bx7eeustzJo1y9PTIfI5LI8kIiIiIqMEQcBf/vIXvPPOO/j444+RmZnp6SkR+STJBm2m9sXyC2gzmmFzZVarcmOyTR0nnd3xUNdRM2LyCaeM500sdfx0dmaMWTYiIpKzuoZmbDpQjdMXb+BGkxohQUokhoXgsfSube2bmprwxBNP4NixY9izZw+ioqI8NGsikmx5pL7qD0eJPzc3dfPIxs/cZFl+7N06gqwjl/IjlkcSyY9c1idbHKmqx/LdZfi69DIAoNnIBtLjhg/A0/fGYUxUH1y8eBGTJ0/G4MGDkZ+fjx49enho5kQESDjTRkRERESOW7enEku2nUaTWgNjX9U3tQdwX5y8hH+X1mF2Uk+seHEG5s6di9/97nfw85Nc3zoinyPZTFvp+xlIWFBi4gzzVPlpiM852OXxsnXazErcTGZWbFG2tv19m2X8favaNBoAEJV13G1zIumTyzfZzLQRyY9c1idraAO2U2hsbbN8cDuhtRmThwJvPz3ZhTMjIltI9qsTcwGbpW0BjAVs5YUpiJt5yK6AzRlt/vXpNgR3puoPRxmUkVpLVZAGVUGa2WPiZh1C3KxDKF2ZYfT5qKzjTg/YKjYko2KDc993IiIiX3Kkqh5Ltp22KWADAEVAILZf7I6j1fUumhkR2crhoE2j0SA1NRUPP/wwAKCiogJjx45FXFwcsrOz0dLS4vAkiYjswfWJiHzZ8t1laFJr7Dq3Sa3Be7vLnDwjIrKXw0Hb3/72N4wY0bHx8csvv4wXXngBZWVl6Nu3L1auXGnXuJc/Hm7yuWEzDonNI2o+GiV2VzTHkT3dYrKPomrTaLEM0JxyM5tMVxSNQUXRGIubgXfOmFmz4XjkoycQ+ahhd0lrMnrxsw8ifnbXzKQx+vvluVrs1KNmNxM3xZoM3fkPknD+A9ObcZPjaj4ahdZhQZ6ehsvWJyIiqatraMbXpZeN3sNmDUEA/nXmMq40NDt3YkRkF4eCturqanz22WdYsGABAO1eHl999RWysrIAADk5OdiyZYtNY+rK9QZMOGPwAbx0dTpKV6cD6Ah+AG0b/M6t8EtXGS/jc4S1JYDmgsPYaUcsBmyANgA7V5wkdqzUD1JtYelauvF15FCSaE2wF/3YMUQ/xi6grhQx+QQCyps8OgdXrE9ERN5i04Fqh8dQANh00PFxiMhxDgVtCxcuxJtvvil2Fbpy5Qr69OkDpVLblDIyMhI1NTWOz5KIyEZcn4jIl52+eMOgrb89mtRtOF1700kzIiJH2N3y/9NPP8XAgQORnp6O3bt323x+Xl4e8vLyAACt6Ei965fqKRTa/60oGoOEaR2bHlvKICXMs7+Mz5r92HTli/rliM7cULt8fSqGTbEts6ZrJhI/+6BYommpJFSp7FjML24Zgdar3WycqWnVH44yeH+snZM97N18XFciyayb/LhqfSIi8hY3mtROGqfVKeMQkWPsbvn/yiuvYO3atVAqlWhqasKNGzfwyCOPYMeOHbh48SKUSiW+//57LF68GDt27DA7ln5LbV2nxphs2+9nMkZ3L5jmthIJj+8XH7+weSQAYPAjJ20ar3JjstPmRiRnnmyp7ar1iYjkwRda/i/ceAhbDl9weJxHUiLwdrbp+/WJyD3sLo9cunQpqqurUVlZiQ0bNuAnP/kJ1q9fj/vuuw+bNm0CAOTn52PixIlOmywRkTW4PhGRr0sMC0Gg0rF+c0FKPySG93LSjIjIEXaXR5qybNkyTJ06Ff/zP/+D1NRUzJ8/36bzLWWxytenWtWUQ2ziYaLMUD/DZqzc0db56cojA7qpfXqTaWdnSvU5e3N0XTMba5rDkDw4uj4REXmLrPRIvL2z1KExBABZaZHOmRAROcQpQdu4ceMwbtw4AMDQoUOxb98+ZwxLROQwrk9E5Iv69wzEvQkD8OWpS3a1/VcogPuGD0C/noHOnxwR2czpmTZHla1tz6bMMp5N6dWz0apxjDUSMdWwwpoMmyW2NCDRb9JhbebQWzg7w6afDbMmw1a64k4AQMKTlj+YW8qwWbOnnq0+rdE21Hk4It2p43pC1abRPp1VJpIbOa1PpPXMuDh8o6pDY6vtG2wHKf3x9Lg4F8yKiOxhdyMSZzJ1o79uzzD9fbfOFqYgMEjbychct0BVfns3xRzrNo52N2u6VHqbs4UpTumeSfIglxv92YiESH7ksj5ZY92eSizZdgqNrda3/w8O8MOrD47AzMwY102MiGwiuUwbERERETmHLvBasu00mtQas6WSCoU2w/bqg4kM2IgkRtKZNh1VfhqUgdr9RvSzbtbQ3z+t9P0MAEDCAvv3cfNlFRuSbX7/7dF5jzfyTnL5JpuZNiL5kcv6ZIuj1fV4b3cZ/nXmMhTQbpytE6T0gwDtPWxPj4tDcmQfj82TiIyTXNBm7/5pUlTz0SijJZzOuo9N101RoYDd4+lvyu0Mciz7JPvI5UMRgzYi+ZHL+mSPKw3N2HSwGqdrb+JGUytCggKQGN4LWWmRbDpCJGEsjyQiIiLyEf16BuKJe4Z5ehpEZCPJZdosUa1JR/ycA2aPqdyY7NQuhqr8NIsNTYw1TXEn/TLQ8vafh7EpiMuUrtZ2V0uYa/7foi+TyzfZzLQRyY9c1ici8h1eEbTp30ulKkiDwk87ZW6y7BhrW9rLqWSV3EcuH4oYtBHJj1zWJyLyHX6engARERERERGZJrlMmzX7q53/QNvsIvox1za7qP5Quxk3Oxlar3y9tjmKnDYMJ/vI5ZtsZtqI5Ecu6xMR+Q7JNSLRlT6aY02w5ox7zKQSrOm6RDqrHFTHFZthyzVYc9d2B0REREREnbE8koiIiIiISMIkl2mLm9WRqSn9vzsAAAmP70flRm3mrE3jZzE7VF6YgmFT5dM5UaFwzbjOzrLJke7fXawTu5ESEREREdlCcve06Vjb2VBqLm4ZgbBJp4w+Z+weOVVBmtUbWzurnLF2ywgAQLiJeUpV1abRAICorOMenglZQy73jPCeNiL5kcv6RES+g+WRREREREREEibZTJsp1R+OMtog5PwHSXZ1k6z5aBQiJkuj4Ygv0t8U3NuVr0+VbSMWe8jlm2xm2ojkRy7rExH5Dq/LtJnq6Ghv+39rArbSFXfaNbY1ytuDFmcf6y2GTj8si4ANkG/nTCIiIiLyLK8L2oiIiIiIiHyJzwRtug257ZHw5D4nzkTbZEVnmA1ZJluOJSIiIiIieZBc0FZemCKWAV7YPBIXNo90yrhNNwNtPkdVkAZVQZpTrq/PG7tiSo3ub1O6MsPxsfKd/zfWp/9vmoiIiIjIVpIL2oiIiIiIiKiDJLpHBsZGIvyPzyA+x/R+ZbqMV/zsgyhblwoAiJtpufGDLotibmySD/1/J+RZcunOxu6RRPIjl/WJiHyHQ5m2+vp6ZGVlITExESNGjMD333+Pq1ev4v7770d8fDzuv/9+XLt2zVlzJSKyGtcnIiIikguHgrbnn38eP//5z3H69GkcOXIEI0aMwBtvvIHx48dDpVJh/PjxeOONNyyOE1h5W8yEnS1MEffuKlubirK12qxa/OyDYvYkbuahLlk2XfYNAKo2jRZ/js856JQsW/n6VJSvT7V8oJc6V2y6UUvp6nSUrk5321z0G7XYSv/fiSm8x8w3OGt9IiIiIvI0u8sjr1+/jpSUFJw9exYKhUJ8fPjw4di9ezfCw8NRW1uLcePG4cyZM2bHChwagcF/egZxsyyXO1qziXbpygwkzC+x7oXoqSga47VNQmq3jAAAhE865eGZEGl5svzImesTyyOJ5IflkUTkbezOtFVUVGDAgAGYO3cuUlNTsWDBAty6dQuXLl1CeHg4ACAsLAyXLl0yen5eXh4yMjKQkZEBzY1b9k6DiKgLZ65PrWh259SJiIiIurA701ZSUoLMzEz85z//wdixY/H8888jJCQEf//731FfXy8e17dvX4v3jfjiN9m6ckSlsg0Rk09YfZ6udDB22hGxjHSol+3fpv8apOxsYYrXvbdS4slvsrk+EZE5zLQRkbdR2ntiZGQkIiMjMXbsWABAVlYW3njjDQwaNAi1tbVi+dHAgQNtGtdUt0dnd4HU3fcWlXXcKeNZQ7+z4ZAp5ks8TZ2nH+h4Y0Chyk9D/DTXdXbU7dtmT3lsZ974/pKWq9YnIiIiIk+wuzwyLCwMUVFR4v0gu3btwsiRIzFhwgTk5+cDAPLz8zFx4kTnzJSIyEpcn4iIiEhOHNqn7fDhw1iwYAFaWlowdOhQrF69Gm1tbZgyZQrOnz+PIUOGoLi4GKGhoWbH0S8/qtyYDACIyT5q77S8ysUtIxDmxuYh54qTbMryEdnL0+VHrlifiEgePL0+ERHZShKba1vzochb7oOyVsUGbXDaersbEuY5XspHJBW6Lwbk8qGIQRuR/MhlfSIi3+HQPm1ERERERETkWnY3InEH/eyaLsPmrr3UygtTMMxCIwrdZtvDZljeX66z2Km+Uf5Jvoflt0RERETOJbmgTb+NfWBQa5fn3VUeaSlgA+wL1oiIiIiIiGzB8kgiIiIiIiIJk0TQ1jI0GOc/SEJ5YQrUzf5QN/sDAAY/chKDHznpkTnVfDTKI9eVo7K1qZ6egkudLUwRM8RERERERM4miaCt29lGRD92DMOmH4bCT4DCz+MNLREx+YRTx3NW4HJh80hc2DzS7vNLV6c7ZR62BClxs+RdRjp0+mFuxE1ERERELiOJoI2IiIiIiIiMk1zQJqj9IKhNT0u1xjmZIneLm3UIpSszULoyw6bzdPu56ThaMpow94Dd5+qTSmapatNoT0+BiIiIiMilJNU9sqJoDJSC2uwx8XOMBx2q/DTE5xw0e+654iQAnmtJnjDf9k20uTWAeVFZxz09BSIiIiIil5Jcpo2IiIiIiIg6SCrTFjvtiMnyR1V+GgCYzKZZyrIB3PTXXvp759mivP08a/a882aObLJORERERGSJpII2AOjWvUX8uXKj9n6umOyjVgVlttDdC8XyOsvsvX9N7sGaDoM1IiIiInIllkcSERERERFJmOSCNrXaH2q1dnPtmOyjiMm2vhFH6ep0q/chi8o6bleWrXM3x85M7V/GLodERERERGQPyQVtRERERERE1EFy97QFBGjsOu9sYQoSptu3B5mqoL3JyWzL981ZasFv6v4vW7N69jb/ICIiIiIieVEIgiB4ehIhilCMVYzv8rh+IxJnkHqXP1Ov19j+cuXrUyX7OogAYK+wCzeEq56ehsNMrU9E5L3ksj4Rke9geSQREREREZGESa48Up9G49yYUoqZqeoPRyHy0RMATGcUje0vJ8XXQkREREREzifpTNuw6Ye77PWlu//M1ao/HOWW6+gCNm9TbqJLJslH2bpUlK1L9fQ0iIiIiHyepIM2IiIiIiIiX+dQ0Pb2229j1KhRGD16NKZNm4ampiZUVFRg7NixiIuLQ3Z2NlpaWpw1VwCA0OqH8vWpYlMRV/HWDJizla01/j53zoB6C0sZVFOv19r9/+QkbuYhxM303jJcT6xPRERERK5gd9BWU1ODd999FyUlJTh+/Dg0Gg02bNiAl19+GS+88ALKysrQt29frFy50q7xS1fcaXi9j0ah5qNRSJhfgmEzDsnmni6plxnGzZLH+6xjKRg39XoT5tq3nQR5hqvXJyIiIiJ3cijTplar0djYCLVajdu3byM8PBxfffUVsrKyAAA5OTnYsmWLUyZKRGQLrk9EREQkF3Z3j4yIiMCLL76I6OhoBAcH44EHHkB6ejr69OkDpVI7bGRkJGpqamwaV5WvbTSSkLPP8HqT3VOuWLtlBAAgfNIpt1zP02WGuvc7PsfyxuJE3sJV6xMRERGRJ9idabt27Rq2bt2KiooKXLhwAbdu3cL27dutPj8vLw8ZGRnIyMhAK5rFx4N6tCCoh+fuMwmfdMptAZunVRSNQXzOQQZsJDuuWp+IiIiIPMHuoG3nzp2IjY3FgAEDEBAQgMmTJ+M///kP6uvroVarAQDV1dWIiIgwen5ubi5KSkpQUlKCAATaOw0ioi64PhEREZGc2B20RUdHY8+ePbh9+zYEQcCuXbswcuRI3Hfffdi0aRMAID8/HxMnTrRp3KaGQDQ1GH5IcmezDlVBmtv2gvO02GlHPD0FIpdw1fpERERE5Al2B21jx45FVlYW0tLSkJSUhLa2NuTm5mLZsmX461//iri4OFy5cgXz58+3adyA4FYEBLfirF6g1r27a8uTzhUn4VxxEgAgfvZBxM/23nLBsxLvRknkDq5an4iIiIg8QSEIguDpSYQoQjFWMR4AULkxGQDQpvHD0PYmHbVbRrj0PjNdwDZkyjGXXcNdzhamiO8bkSftFXbhhnDV09NwmP76RETyIJf1iYh8h93dI10lJvtol8duNQSh9P/uAAD4BavFDX9LV6c7Zf8sOQRrOgzYiIiIiIjkRXJBmz5d1q2bnxrR7YGa7jHAug2Paz4aBcB9WwYQERERERE5k0ObaxMREREREZFrSTrTpl8qef4D7X1nMY91LZ80Rz/DZs1G0rpGHlIpM6z+cBQiH3VelrCiaAy7RhIREREReRHJBW268sfO97ZpNAqHx7ZmE2lHO1Wq1qQjfo7xss2KDdrXFjvV+sDTmQEbwDb/RERERETehuWRREREREREEia5oK21SYnWJiXK16caPB479ahNGSp7hU06hTAbtxfQ3xtNP8um3zQFsP41VBSNsen65Dmlq9ONPq6/95+cVRSN4b9XIiIiIheTXNBGREREREREHSS3ubYl5jaP1mW8NC3+Ju8rczVXbwROQNnaVMTNOuTpaZAFctm8lptrE8mPXNYnIvIdXhe02auiaAxabwUAABIWlLj0WvYqW6ctCdVtHm6Lqk2jAQBRWcedOifyDqWrMpAwT1r/ruXyoYhBG5H8yGV9IiLfwfJIIiIiIiIiCfO6TFv5+lQMm2E8E2VPS31zdA0W2Cbf81gS6X3k8k02M21E8iOX9YmIfIdXBG2q/DSr9lgzprz9PrdhEtks291UBe0bis+27/3zZebunyTL5PKhiEEbkfzIZX0iIt/B8kgiIiIiIiIJ84qgzd4sG6DNsOmybGcLUwz2VHM11Rrje3i5k9DqB6HVOX9mXbMTXyHnLFvZ2lTLBxERERGRJEiuPNLYfWSlq9ORMNczLfxtda44CUOmHLP6+JqPRiFi8gkXzojIM+RSfsTySCL5kcv6RES+wysybURERERERL5KEkFbc0x3qPK1DTNipx3p0q1RP8t2+ePhbp2brWzJsgFglo2IiIiIiMySRNAWWHnb6vvWBkw449rJmHD+gySc/yDJKWPptiYwpXTFnU65jj5VQZrYSdLgWiszULoyw+nXI99TUTRGLG8mIiIiIueRRNBGRERERERExkmuEYkumxX9mG1lhkSuVrtlBAAgfNIpD8/EO8jlRn82IiGSH7msT0TkO5SenkBnpoK1c8XaYK7ldjfEz/GOTpIkLwzWiIiIiMgTLJZHzps3DwMHDsTo0R17dF29ehX3338/4uPjcf/99+PatWsAAEEQ8Ktf/QpxcXFITk7GwYP2769GRGQJ1yciIiLyBRaDtjlz5mD79u0Gj73xxhsYP348VCoVxo8fjzfeeAMA8Pnnn0OlUkGlUiEvLw9PPfWUzRMytQH2kCnHMGTKMauzbLY02DhXnCRm8hx1YfNIp4xDRJa5e30iIiIi8gSLQds999yD0NBQg8e2bt2KnJwcAEBOTg62bNkiPj579mwoFApkZmaivr4etbW1Nk2oTaNAm0aBi+33D9krYX4JEuaXoGJDssVujbqA0BkGP3LSKeMQkWXuXp+IiIiIPMGu7pGXLl1CeHg4ACAsLAyXLl0CANTU1CAqKko8LjIyEjU1NUbHyMvLQ0ZGBjIyMtCKZnumQUTUBdcnIiIikhuHW/4rFAooFAqbz8vNzUVJSQlKSkoQgEDx8bhZhxA36xDCTDR9KFuXatN1YqceRezUoyafV61Jt2k8sl31h6M8PQXyUc5en4iIiIg8wa6gbdCgQWJZUW1tLQYOHAgAiIiIQFVVlXhcdXU1IiIinDDNDnEzDzl1vM73yFlTTtmZrZtuO2uTbnuUG7lf0NUiHz3h9muS7/Lk+kRERETkCnYFbRMmTEB+fj4AID8/HxMnThQfLygogCAI2LNnD3r37i2WKRERuQPXJyIiIpIbi5trT5s2Dbt370ZdXR0GDRqEP/zhD5g0aRKmTJmC8+fPY8iQISguLkZoaCgEQcCzzz6L7du3o3v37li9ejUyMix3cJTz5rXl67XlnMNmGM8Qnv8gSdybrqJoDAAgdtoR90yOyIXcsXkt1ycisgc31yYib2Nxc+2ioiKjj+/atavLYwqFAsuXL3d8VkREVuD6RERERL7AYtAmFWVrtRkrhb+AYdMPe3g21jOVYdPRZdkA12TYVPlpAID4HG4kTERERETkjSQdtOmaZgybfhhxsxxvQKIrVQwMajFojqELCJ1xDalhsOY8uk3fh3rRlwZERERE5P0cbvlPRERERERErmOxEYk76N/or59dM6X0/+4AACQ8vt/1k9Ndc8Wd2ms+uc/isbq93zpvJ+CrytamyjKLSebJ5UZ/NiIhkh+5rE9E5DskVx7Zpu6a/FPlpwGCdoPc+DkHLAZrpavTkTDXuQGTNcGaDoM1QwzYiIiIiIjsx/JIIiIiIiIiCZNUpu1sYQrip3dtnKHfTKNsXSriZprP3Dg7y+ZuFUVj3LJX29nCFK9tqiGlpiD6czlXnAQAGDLlmLlTiIiIiIisJrl72kr/X/v9ak/sR+XGZABATPZRj83NHKl8QK/cmCzZ94jcR2rbO8jlnhHe00YkP3JZn4jId7A8koiIiIiISMIkVR4JAP491QC03RoTsq1v/uEJns6wVWzQZiJjJZ5lc0VjGOpKKhk2IiIiInIuyQVtw2aYv1+tfH2qxWN8RexUaQdrOp4I2Gq3jED4pFNWHy+le+SIiIiIiPSxPJKIiIiIiEjCJJdps5TxGDbjUEdZoJdkmsj9bMmyAcywEREREZF0eU2mrXRVBkpXZQDQBmu+GLCVrk739BSok6pNow1+Ly9MQXn7Fw9ERERERM7gNUEbERERERGRL5JceaSpMjW/AI2bZyI9vtiBUZddlNpr12XYorKOGzw+jGWWRERERORkkgvaTImbxY6RnnC2MMVoIO2uLp6eDNZMTgA+5AAAGK9JREFUvXaga7BGREREROQqLI8kIiIiIiKSMMkFbWcLU8QOkmS78vWpTh3PXBdPuWNHSSIiIiKSAkkEbS1Dg1G5UdvGf+j0w+KHZVVBGlQFaZ6cmtfxhWDKXc4VJ5l8jl0iiYiIiMhdJBG0ERERERERkXEWg7Z58+Zh4MCBGD26Yz+ql156CYmJiUhOTsYjjzyC+vp68bmlS5ciLi4Ow4cPx44dO6yaRLezjYjJPgpVQZpBeWT87IOIn33Q1tdEMnKuOMlsxsuVhkw5ZvK5YdMPs1OkBLhjfSIiIiLyNItB25w5c7B9+3aDx+6//34cP34cR48eRUJCApYuXQoAOHnyJDZs2IATJ05g+/btePrpp6HRsFU/EbkG1yciIiLyBRaDtnvuuQehoaEGjz3wwANQKrW7BWRmZqK6uhoAsHXrVkydOhWBgYGIjY1FXFwc9u3bZ/Vk4mcfNLinTUf/vraKojFWj9eZJ7M2ZJ8hU46ZzXiRb3Pn+kRERETkKQ7f07Zq1Sr84he/AADU1NQgKipKfC4yMhI1NTWOXsKgRLK1SYmKojF2BW++GACwYQb5MnesT0RERESu5tDm2kuWLIFSqcSMGTNsPjcvLw95eXkAgFY0OzINIqIuuD4RERGRXNgdtK1Zswaffvopdu3aBYVCAQCIiIhAVVWVeEx1dTUiIiKMnp+bm4vc3FwAQIgi1Ogxxij8BcROO2LvtK1WsSEZsVOPuvw6rsZmGeSLPLU+EREREbmCXeWR27dvx5tvvomPP/4Y3bt3Fx+fMGECNmzYgObmZlRUVEClUuHOO+902mQBuK2bpCcDtk9rDnjs2uQePzrSIv78cc1+7LhwGDsuMMB2Bk+uT76A65Nv4fpERCQNFjNt06ZNw+7du1FXV4fIyEj84Q9/wNKlS9Hc3Iz7778fgPZm/xUrVmDUqFGYMmUKRo4cCaVSieXLl8Pf39/lL4KIfBPXJyIiIvIFCkEQBE9PIkQRirGK8Z6eBpHL6TJs343pBtWadABA/JwDYldTexvl/LlyD16KyXTOJJ1kr7ALN4Srnp6Gw7g+ETmG6xMRkeMcakRCRLb5bkw3AIBqTTri52jLzKL39gDGOtbV1NkfiCo3JiMm2/vv6SQiz5NawEZE5I0cbvlPRERERERErsNMG5EbfVyzHwAwIaI9wwbg/NhbnpySUZ2zbKp87Qb38TnuaQRERO7XsT7d4eGZEBFRZwzaiNwoUBEAANp72BwsiXQnBmtE8qdbn4iISHpYHklERERERCRhzLQRudHPBqcAAIbAuixb2bpUAEDczEMATDcIUeWnMRtGRA7RrU9ERCQ9DNqIJOrC5pFou6z9v6iqQHtPmaLV+A4d1gZsZWvbg8BZh2yaS+XGZACG97rpd8AkIiIiItdheSQREREREZGEMdNGJDEXNo8EAAx+5CTOFmrLlboFtgIAIh89YfSc0tXp6NuvAQDQolZCoRDEMfQJbQrx5/MfaDf0jn6sa6lm2bpUsSQT6NpNEgCzbERERERuohAEwXi9lRuFKEIxVjHe09Mg8riydalou639LkXZsxVDpx92aLzy9akYNsO2Ukidig3aksjYqUfF4FEZoDEa5BmzV9iFG8JVu64tJVyfiORHLusTEfkOlkcSERERERFJGIM2IgmJm3kIiiANFEEasSRSny7j1ZlqTbrRx4fNOISKDcli1syYs4Up4rg1H41CzUejAGgzbLFTtWWRQ6cfxtDph63OshERERGR8zBoIyIiIiIikjDe00YkIZUbk6Fu9QcAg0YgukyY/j1u1R+OMtmYxB06NyvpTC73jHB9IpIfuaxPROQ72D2SSEKMdWkEYLQhiSsCNlW+dj84a/Z9EzQKi8cQERERkeNYHklERERERCRhDNqIvEj1h6OMPm6qEYmt4nMOWpVl0x1LRERERK7H8kgiL2KqJDIg2LDTpLF74Co2JIvdII35YWsiBk483eVxVX4aAzQiIiKiTuoamrHpQDVOX7yBG01qhAQpkRgWgsfSI9GvZ6BTr8WgjYiIiIiIyEpHquqxfHcZvi69DABoVreJzwUpL+LtnaUYN3wAnr43DmOi+jjlmgzayGuUrsoAACTMK/HwTBz358o9eCkms8vj+lmt0tXpSJh7QPt4e/lj/JwDHceuSRczbDHZRw2ya+qGAPG48vWpAABBDYPHdH1jlQEa7Q9NhnOp3Kjd2y0++6DBef5KjXhNnYqiMYiddsSKV05EUmdqfSIiImDdnkos2XYaTWoNjPXgb2oP4L44eQn/Lq3Dqw8mYmZmjMPXZct/IjNs6abobJc/Ho4BE85YdezZwhSjHSY9SS4ttbk+EcmPXNYnInIvbcB2Co2tbZYPbhcc4IdXHxzhcODGRiRERERERERmHKmqx5Jtp20K2ACgsbUNS7adxtHqeoeubzFomzdvHgYOHIjRo0d3ee6tt96CQqFAXV0dAEAQBPzqV79CXFwckpOTcfAgmxeQd7Olm6Kztaitr172RJatNO8Ot1+zM65PRERE5A7Ld5ehSa2x69wmtQbv7S5z6PoWg7Y5c+Zg+/btXR6vqqrCF198gejoaPGxzz//HCqVCiqVCnl5eXjqqaccmhyRL1MoLFcuny1MwdnCFI8EUAm5+91+zc64PhEREZGr1TU04+vSy0bvYbOGIAD/OnMZVxqa7Z6DxaDtnnvuQWhoaJfHX3jhBbz55ptQKBTiY1u3bsXs2bOhUCiQmZmJ+vp61NbW2j05IiJzuD4RERGRq206UO3wGAoAmw7aP45d3SO3bt2KiIgIjBkzxuDxmpoaREVFib9HRkaipqYG4eHhXcbIy8tDXl4eAKAV9kedRHI1+JGTYufHYTMOGT3Gk81HLO375ilcn4iIiMiZTl+8YdDW3x5N6jacrr1p9/k2B223b9/G66+/ji+++MLuiwJAbm4ucnNzAWi7szlD6WptW3Rdm3Qib6cL1io2aNvvmwqSytenmgzsXEWKAZuU1yciIiLyTjea1JYPsmqcVrvPtTloKy8vR0VFhfgtdnV1NdLS0rBv3z5ERESgqqpKPLa6uhoRERF2T46IyBZcn4iIiMjZQoKcs7V1SFCA5YNMsHkGSUlJ+OGHH8TfY2JiUFJSgv79+2PChAn4xz/+galTp2Lv3r3o3bu30dIjV2GGjeTKVFZLl4ETnPMFkNeT8vpERERE3ikxLASByosOlUgGKf2QGN7L7vMtNiKZNm0a7rrrLpw5cwaRkZFYuXKlyWMffPBBDB06FHFxcXj88cfx3nvv2T0xIl9VtjYVZWtTxY29zYmdehSxU48ibqZrSyMrisagomiM5QPdjOsTERERuVpWeqTDYwgAstLsH8dipq2oqMjs85WVleLPCoUCy5cvt3syRES24PpERERErta/ZyDuTRiAL09dsqvtv0IB3Dd8APr1DLR7Ds4p0CTykLK1qYib5d4GHK5mz+txdSOS2GlHXDY2ERERkdQ9My4O36jq0Nhq+wbbQUp/PD0uzqHrWyyPJJIyuQVs+s5/kGT1sfZu9khERERElo2J6oNXH0xEcIBt4VNwgB9efTARyZF9HLo+gzb6/+3dfUxUZ74H8O/A4GBXaUN1kDvjAsMgAwzDa5G/dG2r25qUVqwsCpGGpvblNv3DRm4TgqUbufTWpLFpTTem20ivVtqQbYdUkG6MWproRcW6DZXOFIYtb0tB5QKNvIw89w8uc2XvDMMMM+cc4fv5azgznvPlOeMv/PKc8xwiIiIiIvKiODcW5duTsCIUENPzL0qiUgErw0JRvj0Jxbmxiz42mzYihfrtru9drzs/TUfnp+keP6sO832qnoiIiIh8U5wbC92Nz5H4mzvQqEMQrp7bToWrQ6BRh+D3yVH4bF9uQBo2gPe0ERERERERLcjFixfReeU8bJ8ew9gUUNfag/b+UYyMTyEiPAym6NV4NlO/qEVH3GHTRnQfMOz5zu32X6ymmRfjEoYhIiIiWqYqKipQUVEBjUYDjQZ4cVO8JMdl00Z0H+j9SwoAQJff5np+W0JJK7RPt0ty/I7/vTQz3kPzSERERLTUXbhwAQ6HAyUlJZIfm/e0ERERERERzUMIgYqKChw8eBBhYWGSH58zbUQy6frMAgCI/cPf3L7/04kMGItnHmmgy29zbU8oaQ1+uH/iaYbN9udsAMCG569IGYeIiIhIUmfPnsXAwACKiopkOT6bNgqan05kAICr8aC5PDVrs4zF1+ConWns4grn/2zXZxav+wsGNmtERES01M3OslVWVkKtlqd94uWRREREREREHjQ2NmJ0dBQFBQWyZeBMGwUNZ9gWz9sMm2tRkj8E95LJ2WfEeVrFkoiIiGgpEkLg4MGDqKysRGhoqGw52LQRycxekznnPrV773Xz1iz5en+b41QaACBu93W377tbJbLz03Q2a0RERLQsWa1WOJ1O5Ofny5qDTRsRERERES1rQ2MTqLvag/Z/jGBk3ImIcDUSo1bj6KH/wL//8Y8ICZH3rjJFN20dJ2cWsogv4mV2tHT982zZvQuKqMPuLmgfHSczFvT/xNMM2yx3q0Rylo2IiIiWquvdwzh6/idcsA0CACac06731CoB56P/hq/++18Q0z2MtPUPyRVT2U0bmzVa7n676/sFfS5UvbDmjoiIiIhmnLjUhaqGdow770KI//++U6iA0DD89cYAmu1DKN9uQnFurOQ5AYU3bURERERERIE207DdwJ2paa+fFQK4M3UXVQ03AECWxo1L/itQd51Z7ggkI/vxLNiPZ7l976cTGbDXZLpWjZx17yWVtmOPuF7PPueNiIiIiGZc7x5GVUP7ghq2e92ZmkZVQzv+1jMcpGSeqYRwNxkorQhVJDaqHpM7BimQt9UOKThsf85e9IOz/0ucxYi4FaBE8mF9Ilp6lkp9IiL/7PvPK/jrjQG3l0R6o1IBv0+Owp+KswMfbB6caSMiIiIiomVhaGwCF2yDfjVswMylkud+HMTNsYnABvOCTRspWtzu65xlW4DZGclAWewsGxEREZES1V3tWfQ+VADqWhe/H194bdpKS0uh1WphNs+9z+r999+HyWRCSkoKysrKXNurq6thNBqRmJiIpqamwCcmmofjVFrAG5j7wXJtbFmfiIiIyBft/xiZs6y/P8ad02jvHw1QooXxunrkc889h1dffRV79+51bTt37hysViuuX78OjUaDX375BQDwww8/oLa2Fm1tbejr68Pjjz8Om82G0NDQ4P0GRLRssT4RERGRL0bGnQHaz1RA9rNQXmfaNm3ahMjIyDnbPvzwQ7zxxhvQaDQAAK1WCwCwWq0oLCyERqNBXFwcjEYjWlpaghDbP8t1Fmah/v55Kv7+earcMRaFl1MuL0upPhEREVHwRYQH5olnEeFhAdnPQvl1T5vNZkNzczM2btyIzZs34/LlywCA3t5erF+/3vU5vV6P3t7ewCQNAKX8Qd/7lxS5I7gVU/A9YgoW9jBnIqW6X+sTERERBZ9pXQQ06sUt6xGuDoEpenWAEi2MX4mdTidu3bqFS5cu4fDhwygoKICvTw44duwYsrOzkZ2djSlIu/oKES1drE9ERETkybNZ+kXvQwB4NnPx+/GFX/ODer0e+fn5UKlUyMnJQUhICIaGhqDT6dDd3e36XE9PD3Q6ndt97Nu3D/v27QMArFmzBrbftGDt2rX+xAmqwcHBgOd64JADtzO9f24+wcgVKErNxly+WWyuFV3yLE4b6Pq0atUq3DY5JMnuq6X63QkW5vKdUrPdr/WJiOS3ZpUGmzesXdRz2rYkrsXDqzSBDzcPv5q2Z555BufOncOWLVtgs9kwOTmJNWvWIC8vD3v27MH+/fvR19cHu92OnJwcr/sbGhpCdnY2rlxR3jLjzOU7pWZjLt8oNZc3ga5PJpNJseOg1HPEXL5Rai5AudmUmouI7g//+jsjmu1DuDN11+d/G64OxSu/MwYh1fy8Nm27d+/G+fPnMTQ0BL1ej7feegulpaUoLS2F2WzGihUrUFNTA5VKhZSUFBQUFCA5ORlqtRpHjx7lymxEFDSsT0REROSrtPUPoXy7CVUNN3BnauHL/68MC0H5dhMs+oeCmM49r03bqVOn3G4/ceKE2+3l5eUoLy9fXCoiogVgfSIiIiJ/FOfGAgCqGtox7rw776WSKtXMDFv5dpPr30kttLKyslKWI7uRlZUldwS3mMt3Ss3GXL5Rai6pKXkclJqNuXyj1FyAcrMpNRcR3T8s+oewKWENbv86ie7bdxAWooJz+v+6t3B1CEJDVHg8SYt3dlqwNXmdbFlVwtdl1YiIiIiIiJaQm2MTqGvtQXv/KEbGpxARHgZT9Go8m6mXfNERd9i0ERERERERKZjsa96eOXMGiYmJMBqNePvtt2XN0t3djS1btiA5ORkpKSl47733AACVlZXQ6XRIT09Heno6GhoaJM8WGxuL1NRUpKenIzs7GwBw69YtbN26FQkJCdi6dStu374taaYff/zRNSbp6emIiIjAkSNHZBuv0tJSaLVamM1m1zZPYySEwGuvvQaj0QiLxYLW1lZJcx04cAAmkwkWiwU7duzA8PAwAKCrqwsrV650jd1LL70kaa75zl11dTWMRiMSExPR1NQUtFxKopT6xNrkG9Ym/3OxNhERKZSQkdPpFAaDQXR0dIiJiQlhsVhEW1ubbHn6+vrE1atXhRBCjIyMiISEBNHW1ibefPNNcfjwYdlyCSFETEyMGBwcnLPtwIEDorq6WgghRHV1tSgrK5MjmhBi5lxGRUWJrq4u2cbrwoUL4urVqyIlJcW1zdMYnT59WjzxxBNienpaXLx4UeTk5Eiaq6mpSUxNTQkhhCgrK3Plcjgccz4XTO5yeTp3bW1twmKxiPHxcdHZ2SkMBoNwOp2S5JSLkuoTa5P/WJt8y8XaRESkTLLOtLW0tMBoNMJgMGDFihUoLCyE1WqVLU90dDQyM2eeer169WokJSWht7dXtjzeWK1WlJSUAABKSkrw5Zdfypbl7NmziI+PR0xMjGwZNm3ahMjIyDnbPI2R1WrF3r17oVKpkJubi+HhYfT390uWa9u2bVCrZxZvzc3NRU9PT1CO7WsuT6xWKwoLC6HRaBAXFwej0YiWlpYgJ5SXkuoTa5P/WJt8y8XaRESkTLI2bb29vVi/fr3rZ71er5g/RLq6unDt2jVs3LgRAPDBBx/AYrGgtLRU8kt9AEClUmHbtm3IysrCsWPHAAADAwOIjo4GAKxbtw4DAwOS55pVW1uL3bt3u36We7xmeRojJX33Pv74Yzz55JOunx0OBzIyMrB582Y0NzdLnsfduVPSeElFqb8za5NvWJv8x9pERKQcst/TpkRjY2PYuXMnjhw5goiICLz88svo6OjAd999h+joaLz++uuSZ/r222/R2tqKxsZGHD16FN98882c91UqFVQqleS5AGBychL19fXYtWsXAChivNyRc4w8qaqqglqtRlFREYCZGZWff/4Z165dw7vvvos9e/ZgZGREsjxKPXc0g7XJN6xN/mNtIiJSFlmbNp1Oh+7ubtfPPT090Ol0MiYCpqamsHPnThQVFSE/Px8AEBUVhdDQUISEhOCFF16Q5dKL2XHRarXYsWMHWlpaEBUV5bpspr+/H1qtVvJcANDY2IjMzExERUUBUMZ4zfI0Rkr47h0/fhxfffUVTp486fqDTaPR4OGHHwYw8wyi+Ph42Gw2yTJ5OndKGC+pKe13Zm3yHWuTf1ibiIiUR9am7ZFHHoHdbofD4cDk5CRqa2uRl5cnWx4hBJ5//nkkJSVh//79ru333k/wxRdfzFnRSgq//vorRkdHXa+//vprmM1m5OXloaamBgBQU1ODp59+WtJcs06dOjXn8iO5x+tensYoLy8Pn3zyCYQQuHTpEh588EHXpUpSOHPmDN555x3U19fjgQcecG0fHBzE3bt3AQCdnZ2w2+0wGAyS5fJ07vLy8lBbW4uJiQk4HA7Y7Xbk5ORIlksOSqpPrE3+YW3yHWsTEZFCybsOysxKWQkJCcJgMIhDhw7JmqW5uVkAEKmpqSItLU2kpaWJ06dPi+LiYmE2m0Vqaqp46qmnRF9fn6S5Ojo6hMViERaLRSQnJ7vGaWhoSDz66KPCaDSKxx57TNy8eVPSXEIIMTY2JiIjI8Xw8LBrm1zjVVhYKNatWyfUarXQ6XTio48+8jhG09PT4pVXXhEGg0GYzWZx+fJlSXPFx8cLvV7v+p69+OKLQggh6urqRHJyskhLSxMZGRmivr5e0lzznbtDhw4Jg8EgNmzYIBoaGoKWS0mUUp9Ym3zH2uRfLtYmIiJl4sO1iYiIiIiIFIwLkRARERERESkYmzYiIiIiIiIFY9NGRERERESkYGzaiIiIiIiIFIxNGxERERERkYKxaSMiIiIiIlIwNm1EREREREQKxqaNiIiIiIhIwf4HnzN9Vp8J9wkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1080x360 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA20AAAEyCAYAAABptTjBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdeVyU9fo//tfAsLkgohLE4ga4giCk1vGUHsvfKUtNzQ0Lc0st006LmZXL0VKPqdXJOpxccCXExBbTT1r2rY674q4sggLiLirKzv37Y5zbAWaYfe57bl7Px8NHMDP3PRcD533mmut9X5dKEAQBREREREREJEsuUgdAREREREREhjFpIyIiIiIikjEmbURERERERDLGpI2IiIiIiEjGmLQRERERERHJGJM2IiIiIiIiGWPSRkRERETkAL169cLXX38tdRgAgNWrV6Nnz55Sh0EmYtJGRERERGSGP/74A4899hiaNGkCX19f/OUvf8GBAwekDosUTC11AEREREREzuL27dt49tln8eWXX2Lo0KEoKyvD77//Dg8PD4fGUVFRAbWab+XrC1baiIiIiIhMlJ6eDgAYMWIEXF1d4eXlhb59+yIyMhKzZ8/GqFGjxMfm5ORApVKhoqJCvC0rKwvdunWDt7c3BgwYgBs3boj3rVmzBi1btkSzZs3wz3/+E61atcLOnTsBALNnz8aQIUMwatQoeHt7Y/Xq1di/fz8effRR+Pj4ICAgAK+99hrKysrE86lUKnz22Wdo06YNmjdvjrfffhtVVVXVfp633noLTZs2RevWrfHTTz/Z5TUj6zFpIyIiIiIyUXh4OFxdXREfH4+ffvoJN2/eNOv4NWvWYOXKlSgoKIBarcbrr78OADh16hQmT56M9evXo6CgALdu3UJ+fn61Y7du3YohQ4agsLAQcXFxcHV1xdKlS3Ht2jXs2bMHu3btwvLly6sds2XLFhw8eBCHDx/G1q1bsXLlSvG+ffv2oV27drh27RreeecdjB07FoIgWPjKkD0xaSMiIiIiMpG3tzf++OMPqFQqjB8/Hi1atED//v1x+fJlk45/8cUX0blzZzRs2BD//Oc/kZycjMrKSqSkpOC5555Dz5494e7ujrlz50KlUlU79tFHH8XAgQPh4uICLy8vxMTEoEePHlCr1WjVqhVeeeUV/Pbbb9WOmT59Onx9fRESEoJp06Zh48aN4n0tW7bE+PHjxSS0oKDA5J+DHItJGxERERGRGTp06IDVq1cjLy8PJ06cwMWLFzFt2jSTjg0ODha/btmyJcrLy3Ht2jVcvHix2n0NGjRAs2bNDB4LaLZqPvvss/D394e3tzfee+89XLt2rc7nu3jxovi9v79/tecDgKKiIpN+DnIsJm1ERERERBZq3749Ro8ejRMnTqBhw4a4d++eeN+lS5dqPT43N1f8+sKFC3Bzc0Pz5s0REBCAvLw88b7i4mJcv3692rE1K2+TJk1C+/btkZGRgdu3b+Ojjz6qtb2x5vM9/PDDlv2gJCkmbUREREREJjpz5gw++eQTMcHKzc3Fxo0b0aNHD0RFReH//b//hwsXLuDWrVv4+OOPax2/bt06nDp1Cvfu3cOHH36IIUOGwNXVFUOGDMH333+P//3vfygrK8Ps2bONXl92584deHt7o1GjRjhz5gy+/PLLWo/517/+hZs3byI3Nxeffvophg0bZpsXghyKSRsRERERkYkaN26Mffv2oXv37mjYsCF69OiBzp0745NPPsFTTz2FYcOGITIyEjExMXj22WdrHf/iiy9i9OjR8Pf3R0lJCT777DMAQKdOnfD5559j+PDhCAgIQKNGjeDn51fnKIHFixdjw4YNaNy4McaPH683IRswYABiYmIQFRWFfv36YezYsbZ7MchhVAJbxBARERERyUpRURF8fHyQkZGB1q1bW3QOlUqFjIwMhIaG2jg6cjRW2oiIiIiIZOD777/HvXv3cPfuXbz11luIiIhAq1atpA6LZIBJGxERERGRDGzduhUPP/wwHn74YWRkZCApKalW8xGqn5i0OanZs2dj1KhRAICcnByoVCpUVFTY5NyrV69Gz549xe8bNWqEc+fO2eTcSqf7eyGqj7g2yRPXJiLn8PXXX6OwsBC3bt3Crl270K5dO6vOJwgCt0YqBJM2ibVq1QpeXl5o1KiR+E93foYt9erVC02bNkVpaalZxxUVFaFNmzZ2iYmI5IlrExERkXwwaZOB77//HkVFReI/e8zPyMnJwe+//w6VSoXvvvvO5ucnIuXh2kRERCQPTNpkavfu3QgKCqp2W6tWrbBz506LzrdmzRr06NEDo0ePRmJiYrX7rl+/jv79+8Pb2xvdunVDVlZWtftVKhUyMzMBAD/++COio6Ph7e2N4OBgzJ49u9pj//jjDzz22GPw8fFBcHAwVq9eDQAoLS3FW2+9hZCQEDz00EOYOHEiiouLq/2sn3zyCfz8/BAQEIBVq1aJ5ywuLsabb76Jli1bokmTJujZsyeKi4vRr18/fP7559WePzIyElu2bKn182u3aa1atQrBwcFo2rQpvvrqKxw4cACRkZHw8fHBa6+9Jj4+KysLf/vb39CsWTM0b94ccXFxKCwsFO9fuHAhAgMD0bhxY7Rr1w67du2q9Zzl5eUYMWIEBg8ejLKyMkO/GiKnwrWJaxMRETkek7Z6Ys2aNYiLi0NcXBx27NiBy5cvi/e9+uqr8PT0REFBAVauXImVK1caPE/Dhg2xZs0aFBYW4scff8SXX36J1NRUAMD58+fx9NNPY8qUKbh69SrS0tIQFRUFAHj33XeRnp6OtLQ0ZGZmIj8/H3PnzhXPe+nSJdy6dQv5+flYsWIFXn31Vdy8eRMA8NZbb+HQoUP43//+hxs3bmDRokVwcXFBfHw81q1bJ57j6NGjyM/PR79+/QzGv2/fPmRkZOCbb77BtGnTMH/+fOzcuRMnT55EcnIyfvvtNwCaPeAzZszAxYsXcfr0aeTm5opvAs+ePYt///vfOHDgAO7cuYMdO3bU6uxUXFyMgQMHwsPDA8nJyXB3dzfht0RU/3Bt0uDaREREdRJIUi1bthQaNmwoNGnSRGjSpIkwYMAAQRAE4ddffxUCAwNrPfbnn38WBEEQZs2aJcTFxQmCIAjZ2dkCAKG8vFzvc/z++++CWq0Wrl69KgiCILRr105YsmSJIAiCUFFRIajVauH06dPi42fMmCH85S9/Eb8HIGRkZOg999SpU4Vp06YJgiAIH330kTBw4MBaj6mqqhIaNGggZGZmirf973//E1q1aiX+rJ6entXib9GihbBnzx6hsrJS8PT0FNLS0mqdt7i4WPDx8RHS09MFQRCEN998U5g0aZLeOLWvUV5ennibr6+vkJSUJH4/aNAgYenSpXqP37JlixAVFSUIgiBkZGQILVq0EH7++WehrKys2uNmzZolPPfcc8Ljjz8uTJkyRaiqqtJ7PiK549rEtYmIiOSDlTYZSE1NRWFhIQoLC8VPhm0pMTERffv2RfPmzQEAI0eOFLchXb16FRUVFQgODhYf37JlS4Pn2rdvH3r37o0WLVqgSZMm+Oqrr3Dt2jUAQG5uLtq2bVvrmKtXr+LevXuIiYmBj48PfHx88Pe//x1Xr14VH9OsWTOo1Wrx+wYNGqCoqAjXrl1DSUmJ3vN6enpi2LBhWLduHaqqqrBx40a8+OKLdb4WDz30kPi1l5dXre+LiooAAJcvX8bw4cMRGBgIb29vjBo1Svw5Q0NDsWzZMsyePRt+fn4YPnx4tQYNe/fuxbFjx/Duu++yTS85Na5NXJuIiEgemLTJVMOGDXHv3j3x+8rKympvJExVXFwsbq3x9/eHv78/li5diqNHj+Lo0aNo0aIF1Go1cnNzxWMuXLhg8HwjR45E//79kZubi1u3bmHixIkQBAEAEBwcXOuaEwBo3rw5vLy8cPLkSfEN4K1bt8Q3IXVp3rw5PD099Z4XAOLj47F+/Xrs2rULDRo0wKOPPmr0nKZ47733oFKpcPz4cdy+fRvr1q0Tf05A8zr88ccfOH/+PFQqFaZPny7e17dvX8yYMQN9+vSpttWLSAm4Nj04lmsTERE5CpM2mQoPD0dJSQl+/PFHlJeXY968eWa3wwY0n5S7urri1KlTSEtLQ1paGk6fPo2//vWvWLNmDVxdXTFo0CDMnj0b9+7dw6lTp2o1A9B1584d+Pr6wtPTE/v378eGDRvE++Li4rBz504kJyejoqIC169fR1paGlxcXDB+/Hi88cYbuHLlCgAgPz8fO3bsMBq/i4sLxowZg3/84x+4ePEiKisrsWfPHvG1ePTRR+Hi4oI333zT6CfZ5rhz5w4aNWqEJk2aID8/H//617/E+86ePYtffvkFpaWl8PT0hJeXF1xcqv9P6Z133sHIkSPRp08f8VNwIiXg2qTBtYmIiByJSZtMNWnSBMuXL8e4ceMQGBiIhg0b1urYZorExES8/PLLCAkJET/N9vf3x2uvvYb169ejoqIC//73v1FUVAR/f3+MHj0aL7/8ssHzLV++HB9++CEaN26MuXPnYujQoeJ9ISEh2LZtGz755BP4+voiKioKR48eBaDpaBYaGooePXrA29sbTz75JM6ePWvSz7B48WJERETgkUcega+vL6ZPn46qqirx/pdeegnHjx+36eDYWbNm4fDhw2jSpAn69euHQYMGifeVlpbi3XffRfPmzeHv748rV67g448/rnWODz74AAMHDsSTTz6JGzdu2Cw2IilxbXqAaxMRETmKStDdV0FUQ1VVFVxdXXH+/HmEhIRIHY5ea9asQUJCAv744w+pQyEiB+HaRERE9Yna+EOoPjtx4gQ8PT3h7+8vdSh63bt3D8uXL8fkyZOlDoWIHIhrExERycG1olKkHMrDmUu3cbukAt6earT398YLMUFo1sjDZs/DpI0M2rx5MyZMmICFCxfKcpbPjh07MGjQIDz55JMYOXKk1OEQkYNwbSIiIqkdzS3EF7sz8Vu6phlXacWD7fGe6ktYujMdvdq1wOQnQtEl2Mfq57Pb9sjt27dj6tSpqKysxLhx4/Duu+/a42mIiMzCtYmIiIissW5vDuZvO4OSikrUlUmpVICn2hUzn2mPUT1aWfWcdknaKisrER4ejp9//hlBQUF45JFHsHHjRnTs2NHWT0VEZDKuTURERGQNTcJ2GsXlVcYffJ+XmwtmPtPBqsTNLt0j9+/fj9DQULRp0wbu7u4YPnw4tm7dao+nIiIyGdcmIiIistTR3ELM33bGrIQNAIrLqzB/2xkcyyu0+Lntck1bfn4+goODxe+DgoKwb9++ao9JSEhAQkICACDtUBoawLva/WVtvKp9736uWHN7W88Ht2WVVHtMaasGAACPnHuwVPn981dVqeCRXSye15pzyoH29XQ/V4zS1pqvtT9fXcrbesKtxussV6WtvUz6mexJ93Wu70pwF2WC+fO77MmUtQkwvj4RkXOT4/pERPL3xe5MlFRUWnRsSUUllu/OxFejYi06XrJGJBMmTMCECRMAAN4qX3RX9al2f/ZHXeDlVQYAuHPbC2EvHQYAZM6NRuioI5qv1z/4OntjF7iqNVlvyAvHTY4j4/PuCJui86btnM6dqvv/Pf/g69yUzggecsLk88tG9v3/qgDk6HxtRO6i6j9vbkpnAJDna5ADk34mU13YFGHW3xIA8XXO/iYSANB6+DHbBeRk9gm7pA7BYsbWJyJybs68PhGRNK4VleK39Kt1XsNWF0EAfj17FdeLSi3qKmmXpC0wMBC5ubni93l5eQgMDDTrHCoXAcUlbgAAzwZlyFjTFQAQNuowbv4YBgCouuwqPr71iKPIuf9GOXPdg2TuwqYIAEBZqRq+PncBAL7PpouPDRv2IGHLWh8t/iK0x9cle2MXAEBFmSvC4g+Lz2fKG/0rW9sDAPwGnEH6qhgAQPjLhzTnTYqs9mb/3IYoAECbkWnibRe3dETp/den9YijhmNMqp48pH/VDeET94v3a5OwykoXtBr24Dm1r4+7a4X4/JXlLmjoWiL+nNrXqux+HGr3CqNJSt7mTigu0vyhunlWwNVVk2jXTALzNncCAJQUazrD1fX70Pf66CbXmeuiAe3v9cW6f695mzuhtFTz87R9IQ3p/3kEAKDyqhQ/OEj/shug1pwwfPwB8djMtdHi+SvuvyYkL7ZYm4iIiKj+STmUZ/U5VABSDufhlcfbmn2sXa5pe+SRR5CRkYHs7GyUlZUhKSkJ/fv3t8dTERGZjGsTERERWeLMpdvV2vpboqSiCmcK7lh0rN1a/m/btg3Tpk1DZWUlxowZg5kzZxp8rGfbQIQsfAVVlS5QuWjCaTXsGDLXRgMAXNVVULtp9o8W3/ZE+JiD4rG6FSvdqpK2UhM0+KT42Kz7FZm2OhUZQFNhA4CqcheEjdZUuzLWdBUrK9rvASDspcNI/6obAIgVq+yNXapVu/RVfnSrMOkrYhE+9sHPYAvan6FtnPEKoRR0Xz9dhn4nxo4zV82Ko7lq/j2QcfuEXbgt3JA6jFrMWZsAbo8kUiK5rk9EJF9jEg/glzNXrD5Pn/Z+WBH/iNnH2S1pM4dn20AEfTQJXl5l4pbIynJXMcnJ2hAFFxdNZlte7IZGPpomD5WVLlCpBPHrlkM12xIz10ZD7a5J8sruac4X/vIhcStgeZm6WkJVkNoBABAw8LTe+DLWdIVXQ80Fy6Ul7kYTI1MSBO22zfIytd7zabcn6m5ZNOR8coT4s9dUM5nL+SbS6DnzNndCeblm62ldP4O+GLXPp/ucpspcpznWlK2ptmbouc35PQDApft/S/4G/pacXfqX3RA+ab/xB0I5b4qYtBEpj1LWJyJynGnfHEFq2kWrz/N8VCCWDosy+zi7bI8kIiIiIiJSivb+3vBQW5c6eapd0D6gsUXHyqrSVlHuCs8Gmo6RFeWuqBI0rQDbjkwTm37U1XRDK2t9dK0qj77batJW4gCgtMRNrLpYszXu4hbN0N6Hnz9lcqy6ajY2MbfyI2fa7a8AoHIBKos11b3wcca3juZ/q9n+GjjopFl/G6YwtmXTHLpNceobpXySzUobkfIoZX0iIse5VlSKvyz8xarr2jzULvjf9L/Jp3ukudzPFYtvuLXXMHk2KENFhSab1b1mTDfJyVwXDZf718Dpbnd0cakS33hDm/jVSIx0r4XTMtTG3ljClrUhyuAbfEHQ34P+nJgY1P2GvqTIvdr3+pI1UzpWZqzWdKjUXrOnpZtUGtsmqqvmdXzG6Es2dTs5XtgUIf6Ozm2IgrtHOYDq1yTqChz04HZDcZg6niBjdUy116WqzLWORz84BgA8GpbV+drX14SNiIiISEmaN/LAE+Et8PPpyxa1/VepgN7tWliUsAHcHklERERERGTUq71C4ak2/uG+Pp5qV0zuFWrxc8ui0lbW1hO5izqjtMQNYaPuz8JaGQvcbzIS/vKhag01tNvhQkccEStFwIOqUeXdB10oDW091FbYzidHoPS2JuPV3ZanraQAtatTNdXqRnm/iubuXgEYqKC20VOZS1+hmZDu6lUhNl4JH/7gudNXxsLLWzMnTbd6pG3GAhiuuun+DLqNUnS3bWorbDXnxGmd2xAlxm3SNtX7r4NQpULosLorTpUVLuLvtc2I2q+N7ty1vM2d9FbgDMWt6+KWjrW2qqrU1X9J2tdKt+Mn8KB5TMgLx43+TYgxbezyoPPpTa9qc92IiIiIyHl0CfbBzGfaY96Pp1FixjZJLzcXzHymPSKDfCx+bnlc0xb6MFotegXBQ06Ig7Ob9svQez2Y7rY83U6Il1I76O3Yp92WV1nhKiY3+hImY7SJTlWla61EsOabe0PdI3W3CBq7bsrQdVraLodV5S7iMG5T6BtD4CjXvg9H8+fSHf68ltBNDg0x5frImo8H5DuOwV6Ucs0Ir2kjUh6lrE9EJI2/vzYf6Y06Q3BR17lVUqXSVNhmPtMeo3q0suo5ZVFpIyIiIiIikrtdu3bh1Hf/wTc/70HigQL8evYqVEC1ypun2gUCNNewTe4ValWFTUsWlTbtJ9mZ66LFJhDhYw5WayShO0dLX1ML3a175Hzyv+1UrbkIOT+lfJLNShuR8ihlfSIix7p79y4iIyPx2WefoV+/fgCA60WlSDmchzMFd3C7pBzenm5oH9AYQ7oGWdx0RB9ZJW26rmxtj7IKTQIXNPikuL3P06vMYPt8U2Ws6QoX17qveQPuX1cHTQKpj24iaWqnQoAJCtUPSnlTxKSNSHmUsj4RkWO9+eabuHTpEtavX+/w5+b2SCIiIiIiojrs378f69evx/HjdY/ZshdZVdpqzv46n6zp1tdyaPUX59L9eWL37nmYtSXSUIMQXfqe81JqB5SWqfXGUhfdRimOUpDaoc45a+bOV6uprqYapnRvpPpDKZ9ks9JGpDxKWZ+IyDHKysoQExOD9957DyNGjJAkBlkkbR6tg+A/awo8G5eiqkozjFoQ9CdX+lq2WyMjsSvUHhUAqj+foS2M5zZEoVFDTdt93cHc9Zmh5NqW0v/7iNguv2a3TqqbOVt3bUkpb4qYtBEpj1LWJyJyjLlz5+LAgQP47rvvoFKpJImB2yOJiIiIiIj0OHnyJD7//HMcOXJEsoQNkEnS5pFzT2z2odvcI3Ot5muVC8Rh05V3XWzynNotlmEDD4u3pa+KEWefBQ46qXdWmu52zLzNnQBA76DnmmoOvTZlq2ZNhmaIZW2IMjjvrabspEhUVWleQ1OPMcaeFTYtr6bF4tf2rrLp607qzBxdYSMiIiJSgsrKSowdOxbz5s1DUFCQpLHIYnukvu1HeZs7oeJ+90jdYdSA7ZINuUlfFQMACH/5kElJnZQDs51JzYSZHEMp24+4PZJIeZSyPhGRfS1btgypqan45Zdf4OJim8KRpWRRaSMiIiIiIpKLc+fOYd68edizZ4/kCRsgw6TtwiZNU4uQwcfFyhMAQNDsIW0bp3+4tiPZ+vm12zDDRxwSb9NW2AxtiQTMq7CJr6uMK07mVA7zv9VsTTVl3l3IC8eRsVrztxQ2+pCRRxMRERFRfSYIAiZMmIDp06cjLCxM6nAAANKnjURERERERDKxatUqFBYW4o033pA6FJHsKm1lpQ9C0jYFyVofXW0uWOX9a92kYmmFzVBr/PJbHgaPsaaJhLapCwCEviD/FvnGKmza6xrbjkwzqcKmixU2IiIiIjLm4sWLePfdd/Hzzz9DrZZPqiS7Spuvz134+txF1oYo5KZ0Rm5K51qDnFUqASpV9f4p2sYdptB9rLYDZF20WwsNSV8RK36tmyjpylwbjZZDj+vttBg+cT/CJ+5/cL6VD86n3YppidBRR8R/tpKdFFnna53zTaRVMQOo1nRGl1ClglBVvdWqdgaZtbI3dhG3qdqSKX9fSmPt75+IiIhICoIg4NVXX8Urr7yCLl1s/77QGvJJH4mIiIiIiCSyefNmnDlzBklJSVKHUovFLf9zc3Px0ksv4fLly1CpVJgwYQKmTp2KGzduYNiwYcjJyUGrVq2QnJyMpk2b1nkuQy21TWlMkbGmKwAg7KXDBh+jZayN/rkNUWyfX09p5/b5DzwtcSTKIWVLbUesT0TkvNjyn4hqunHjBjp37oyUlBQ89thjUodTi8VJW0FBAQoKCtC1a1fcuXMHMTExSE1NxerVq+Hr64t3330XCxYswM2bN7Fw4cI6z6X7pshQZ8YrW9sDAPwGnKl2e9Z6zXZE3S2Ul1I7mPzm25JEzZxE0RQXt3QEADz8/CmbnM8c5gwIr0l3EDpRTVK+KbLX+kREysCkjYhqGj16NJo0aYJPP/1U6lD0sviatoCAAHTtqkleGjdujA4dOiA/Px9bt25FfHw8ACA+Ph6pqam2iZSIyERcn4iIiMhUO3bswG+//Yb58+dLHYpBFlfadOXk5ODxxx/HiRMnEBISgsLCQgCai/maNm0qfm+I9pPsjM+7I2zKvjofez45Qm8zDwDISLxfAYvXXwHTrSqZMxPMmWRv7ILWI45KHYZs8fVxHLl8km2r9YmIlEMu6xMRSe/OnTuIiIhAQkIC+vbtK3U4BlndiKSoqAiDBw/GsmXL4O3tXe0+lUoFlUql97iEhAQkJCQAAMpRCgDVEras9dFo0LAEABCgs9Wx9Lbh9vhqj4o6Y9XdAmhKsmZsG6SttzVquxfWlVRc/a4dAKBF/7N675cyIclcF232VklzhmRbKmN1jNjy39DrI/XAdrIPW65PREREpDwzZ85E7969ZZ2wAVYmbeXl5Rg8eDDi4uIwaNAgAMBDDz2EgoICBAQEoKCgAH5+fnqPnTBhAiZMmABA80k2EZEtcX0iIiKiuvz5559ISUnBiROWz0V2FIu3RwqCgPj4ePj6+mLZsmXi7W+//TaaNWsmXuh/48YNLFq0qM5z6W4/0jYWqSp3EasjuSmdzRoynb4qRhzMrU997xKpnTsX8oL+baZEtiDl9iN7rU9EpAzcHklEJSUliI6Oxrx58zB48GCpwzHK4qTtjz/+wF//+ldERETAxUXTz+Sjjz5C9+7dMXToUFy4cAEtW7ZEcnIyfH3r/qTaK/RhtF48AYGDTurtSGhK0lbXtW5ylpHY1eA1eEqUnRQpjlzI29xJb9fKrPXRcHGtBGB4PIM1Mtfe/xt7kV0v7UnKN0W2XJ+YtBEpD5M2Inr//fdx+vRpbN68WepQTGLx9siePXvCUL63a9cuiwMiIrIW1yciIiIyJC0tDQkJCTh61Hma01nc8t+WqqpUKClxA6CpsIWOOiI2AQGA0vv3AZqmEvq0HHoc+d92EhtbSEW7vdNUcqyyabdPAkD6yliLz5OdFCkONNfSrZwZmg3XNu4IWg8/ZpcqG6CpsElRZUtfpf9vl4iIiIgco6KiAmPHjsXChQsREBAgdTgms0nLf2uZs/1ItxMg8CCpCB9zULzN3Lbuth6W7QimXpenxOvX2LbfOShl+xG3RxIpj1LWJyIy36JFi/Dzzz/j//7v/wx2kZYjq1v+ExERERERyV16ejoWLVqEAwcOOFXCBsik0ubROggBc1+ttVVQd4uksSrYpdQO8NeZ56aPo6tOmWujq23DO5+seX5LG6U9zo4AACAASURBVKaY2kVTu0WzbRwbbZB0lPJJNittRMqjlPWJiExXVVWF3r174/nnn8e0adOkDsdssqi0eeTc05uweTXUDLWtee2T9jop3WueSsuM/yiO3iIY+uKRakObtcla1oYotLVg5ICpYw/qStZyvolU9ABpUwaU16W+j4MgIiIiUqKEhASUlpZiypQpUodiEVkkbURERERERPaQm5uL999/H7/99htcXV2lDsciskraclM6i1+HDTlssBNjVWXtF7vl0OM4tyEKACStlOjOIQOgt6pVs8qWdT9uS6pv5rJHlU3b0VO3QYxUrG1QwiobERERkXIIgoBJkybh9ddfR6dO0naZt4askjbd7X/pX3VDeNx+vY8ztP2vUcOSWrfpbk/UZa+Okea2qb/6XTu07e/ciYIckjUiIiIiopo2btyI8+fP49tvv5U6FKvIKmkjIiIiIiKyhatXr+If//gHvvvuO7i7u0sdjlVk0T1StzubuY0kMtdptlCGjnpQfcvb3Mng4GZ9Lm7pCAB4+PlTJh+jZGzG8YB2y66pTWDoAaV0Z2P3SCLlUcr6RER1i4uLQ0BAABYvXix1KFZjpY2IiIiIiBTlhx9+wL59+3DsmDK6psuu0paRqLnWTO1eWa3aptusI3Pt/erai+bNIUtfEQsACB970OqYyTGsnW2njyMbv9RnSvkkm5U2IuVRyvpERPrdvn0bnTt3RmJiInr37i11ODYhu6TNXM6+fS1zXbS4tdPYz5K1PpoDs8lpKOVNEZM2IuVRyvpERPpNmjQJlZWVSEhIkDoUm+H2SCIiIiIiUoTffvsN33//PU6ccM6CjiEuUgdQ04VNEbiwKcLkxwcPOWFylS1zXbTYuMQYQzPispMiTY7NEO2WP6B6A5WqKhWqqlQGj6urypa1IUrc9udoeZs7IW+zbedeZG/sIjalAR6MbtCnILUDClI7GLz/nESvCxERERE5TnFxMcaNG4fly5fDx8dH6nBsSrbbI89tiILKpQpA9dlnWRui4KLndjnK3tjF6mHP5FhZG6J4rZuNKGX7EbdHEimPUtYnIqpu+vTpOH/+PJKSkqQOxea4PZKIiIiIiJzawYMHsXr1ahw/brvmdXIi26TN0Jwwd/cKvdshc76JRKth9qm8Za6NNqtTpbmz5hx9Pmei3XZpaO6erauZrLIREREROZfy8nKMHTsWn3zyCfz8/KQOxy5kuz1S18UtHSEImmu9qqpUknSK1F5TZa/EkBxPe60bB4nbh1K2H3F7JJHyKGV9IiKN+fPn488//8SPP/4IlcpwfwhnJttKGxERERERUV1Onz6NZcuW4dChQ4pN2AAZJ226WxIffv6UWcde2BSBkBdq72e1dCg3UL8rbNqB52Hxhy26X65YYSMiIiJyXpWVlRg7dizmzJmDkJAQqcOxK9kmbXUlVsa2KupL2M4nRyB0qGWDqc8nR6DlUNtd1HhxS0ezE1FjtKMMdEcImMKULYLaZCw7KVJvx057JGsZa+4ngi85VyJIRERERI6xfPlyuLq6YuLEiVKHYndWz2mrrKxEdHQ0nn32WQBAdnY2unfvjtDQUAwbNgxlZWVWB0lEZAmuT0RERMqUk5ODOXPm4Ouvv4aLi+xGT9uc1T/hp59+ig4dHgw2nj59Ot544w1kZmaiadOmWLFihUXnTV8Ra/C+VsOOiVW2/G87If9b44OdramUtRx63OTh1bqDs2u6uKWjSVW23JTO1b7PXBstbu00JHTUkVpVths/hBuJVlNhM3WboCPn4oW9dNiiKltuSudar19N55Mj6vw9kfUy1nRFaesGUodht/WJiIiIpCMIAl555RW89dZbaNeundThOIRV3SPz8vIQHx+PmTNnYsmSJfj+++/RokULXLp0CWq1Gnv27MHs2bOxY8eOOs+j253tytb2AAC/AWdwYZPmjXXIC8erbePTJk+G2rMb2sbnTKy5/s5UNa/9u7ilIwDzryEk0kfq7mz2WJ+ISBmkXp+IyDqJiYn49NNPsW/fPri5uUkdjkNYdU3btGnTsGjRIty5cwcAcP36dfj4+ECt1pw2KCgI+fn51kdJRGQmrk9ERETKc+nSJbz99tvYsWNHvUnYACuSth9++AF+fn6IiYnB7t27zT4+ISEBCQkJAIBylIq3+w04I35dXqYJL2t9NNqOfFBxMjYA2ZoqmymdEK9+pynDtuh/VrzN0kYg+pzbEIXQkeadJztJ05yl9fBjJlfMykoe/PqvbG2Pe7c9zIzUsJrNW4xVR62h3Q5p7vy+S6mabXP+A0/bPCaSlr3WJyIiIpLWlClTMG7cOERH133pkNJYvD1yxowZWLt2LdRqNUpKSnD79m08//zz2LFjh1Xbj9JXxQAAwl8+ZElYtWRv7AIAKL/lgfCJ+8XbLd0KmLUhyi6JB5HSSLn9yF7rExEpA7dHEjmnb7/9FjNmzMDRo0fh6ekpdTgOZXEjko8//hh5eXnIyclBUlIS/va3v2H9+vXo3bs3UlJSAGj2mw4YMMBmwRIRmYLrExERkbLcvHkTU6ZMwYoVK+pdwgbYYU7bwoULMXz4cLz//vuIjo7G2LFjzTreWIUtfUUswsceNHoesWI3Qv/5dCts2mpc6xFHjZ7XUJVNu0Wvqkpl05luzkZ3m6at6TapsQVbbmkl52Dt+kRERETSeOutt/D888+jZ8+eUociCZskbb169UKvXr0AAG3atMH+/fvrPoCIyEG4PhERETm3nTt3YufOnThxwrz+BUpiVct/W9G9ZsRYpcaaa8rSV2pmv4WPMV6ps6fclM5i04xzG6JMnpNWH5lTBQWAjNWaCmvYaOuviazZTMUWdlzU/K7/v4eNz/yTu4LUDgioo4mLUq4Z4TVtVF8oaX0yRinrE1F9cPfuXUREROCLL77A008/LXU4kpFd0qYr5xtNAqcdpA3cT+oEFYC638jbc5ueLThiDpujmTI4nOoPpbwpYtJGpDxKWZ+I6oM33ngD169fx5o1a6QORVI2v6aNiIiIiIjIWnv37kVSUlK93hapJeukTVthu7ApAiVF7gCA8OGmbXvTVth0Z3jpq9xJxRkrbJnrouts2mGrKlt2UqRsK6REREREZH+lpaUYO3Ysli1bhmbNmkkdjuRktz3y3P0hzEq4zit7Yxe9Wzgz10bbJGnTvlaNGpZY3FHRlteAAcrc9kmWUcr2I26PJFIepaxPREo2a9YspKWlITU1FSqVSupwJCfrShsREREREdUvx48fx/Lly5GWlsaE7T7ZVdqMSV8Za7T7o253Rlsw5XxZ6zUVprZx0lSYdLeBXvs+HADQ/Ll0SWKpD6T+fTsDpXySzUobkfIoZX0iUqLKyko8+uijGD9+PMaPHy91OLLhFElbzjeR4nVouSmdUVriBsB2Q5GvftcOANCi/1mbnM9Z6L6udeGWR7KEUt4UMWkjUh6lrE9ESrRkyRL88MMP2LVrF6tsOrg9koiIiIiIJJeVlYWPPvoIe/fuZcJWg+yStotbOgKo3olQtxoUPOQE0lfF2PQ5DVXYlNQURR9Tu2iaU2HjtkEiIiIiMpcgCBg/fjxmzJiB0NBQqcORHdklbdqtj3UJf9l4p0NbJFxySdbslTxmrY+2eXKl1GTN1K2kRERERGS+FStWoKioCFOnTpU6FFmSXdJGRERERET1R35+PmbMmIFffvkFajXTE31k96rozjU7nxwBAGg59DgubNJ8rVIJxjs5bohCW5lUyWyhUcMSu5xXqVUxW5LTQHYiIiIipREEAZMnT8bkyZMREREhdTiyJbukTStrfTTaDn2QVIS8cNzkY6VM2OpKGDPX3e/CqNP1MjspEq2H150QaAdn22qUgTYBNuc1lYP0lbEAYHTkgy0xWSMiIiKyn02bNiEzMxPJyclShyJrsk3aiIiIiIhIua5fv46pU6diy5Yt8PDwkDocWXOKOW26LmyK0FshKkjtgICBp81+7uyNXaptySTHUlK3SaVty7WWUuYgcU4bkfIoZX0icnYvvfQSmjVrhqVLl0odiuy5SB2AuQxt6bMkYQNgUsKmTSzsQTviwBTXvg+3WxxSaRt3RBEJGyDttlwiIiIiZ/LTTz/hjz/+wLx586QOxSlweyQRERERETnMnTt3MHHiRKxYsQINGzaUOhyn4HSVNktpm4BYwtaVoBs/PKiY6Q4RN6b5c+k2jYOIiIiIyNFmzJiBJ598Ek8++aTUoTgN2VXaspM0LdZbDz8m+fVO9uqy6Pssky9raX83gqAZCWGN88kRVp+jLrp/00RERET12e+//44tW7bgxAnrO6LXJ/Wm0kZERERERNIpKSnBuHHj8O9//xtNmzaVOhynIovukZ5tAxG8YGKdM7EyVscAAMJGH8K5DVEAgDYmNH5w1plkZJnclM4AYJN5dmQdpXRnY/dIIuVRyvpE5Gzee+89ZGRkYNOmTVKH4nSsqrQVFhZiyJAhaN++PTp06IA9e/bgxo0beOqppxAWFoannnoKN2/etFWsREQm4/pEREQkH0eOHMGKFSvw+eefSx2KU7IqaZs6dSr+/ve/48yZMzh69Cg6dOiABQsWoE+fPsjIyECfPn2wYMECo+dxP1csVtnSv+qG9K+6AQByvolEzjea64HCRh9C2OhDADQVtppVtuyNXR58ff8aIkBTYbNFlU03FiXSVjL1yVofbdexBzVpq6OWCB5ywmiVLXNdtFWNacg52Gp9IiIiIuuUl5djzJgxWLRoEfz9/aUOxylZvD3y1q1biIqKwrlz56BSqcTb27Vrh927dyMgIAAFBQXo1asXzp49W+e5PNoEIuijySZtd8xOijTa0MHSIceGBnc7A24DJbmRcvuRLdcnbo8kUh5ujyRyrAULFuDXX3/F9u3bq/3/MpnO4kpbdnY2WrRogZdffhnR0dEYN24c7t69i8uXLyMgIAAA4O/vj8uXL+s9PiEhAbGxsYiNjUXVnbuWhkFEVIst16dylDoydCIiIkU5e/YsFi9ejP/85z9M2KxgcaXt4MGD6NGjB/788090794dU6dOhbe3Nz7//HMUFhaKj2vatKnR60bq4yfZ2spYWYkaoS+aPtJAu0Wz1TDpRyJYylla4JtS1SXDpPwkm+sTEdWFlTYix6iqqsITTzyBoUOHYsqUKVKH49QsntMWFBSEoKAgdO/eHQAwZMgQLFiwAA899BAKCgrE7Ud+fn5mnddQ97+LWzoCMG8YdV3yNncCAAQNPmmT85lCt+ulOdsYdTtn6nbYdLZkDdC87q0H2y8RSl8ZCwAIH3PQ6nMxYXNe9lqfiIiIyHRfffUVqqqqMHnyZKlDcXoWb4/09/dHcHCweD3Irl270LFjR/Tv3x+JiYkAgMTERAwYMMA2kRIRmYjrExERkbQuXLiAWbNm4euvv4arq6vU4Tg9q+a0paWlYdy4cSgrK0ObNm2watUqVFVVYejQobhw4QJatmyJ5ORk+Pr61nke3e1Hutv/6oMrW9vDb8AZhz1fxuoYsQsnkT1Jvf3IHusTESmD1OsTkdIJgoB+/frhL3/5C2bOnCl1OIogi+HaprwpKkjtAAAIGHjaESHZnXa7570iD7OuaSOSu5xvItFq2DHFvCli0kakPEpZn4jkat26dfjXv/6FgwcPws3NTepwFMHia9qIiIiIiIh0XblyBW+++SZ+/PFHJmw2JOukTbdZiLbC5qiOfqbMbNNtLGIuWzVUIZKb+rK1mYiIiGp7/fXXER8fj9jYWKlDURTZJW2617SVl9e+aNFRHf1M6e5oSbJGRERERKREW7duxaFDh7By5UqpQ1Ec2SVtRERERETkXG7duoVXX30V69atQ4MGDaQOR3EsbvlvS+VtPVGQ2gG5KZ3h7l4Bd/cKAJqqmlSzsrTzvsh62mHaSpWR2BUZiV2lDoOIiIhIMu+88w6effZZ9OrVS+pQFEkWlTa3rBLxmjXtdWJSs8VwZl0ZiV0RFn/Y6vNkro0GAIs7Tmatj7bJUG5zri1U+pBqW/xeiYiI7O1aUSlSDuXhzKXbuF1SAW9PNdr7e+OFmCA0a+QhdXjkxH799Vds27YNJ06ckDoUxZJF0kZERERE9nE0txBf7M7Eb+lXAQClFVXifZ7qS1i6Mx292rXA5CdC0SXYR6owyUndu3cP48ePx/Lly9GkSROpw1EsWWyP1FVZ7oLKcsNhyaUSZ66w+MPIToo0e6tg5rroat+HvnjEqrlutqiyAfKpnjnr3wMREZEjrNubg+H/3YufT19GaUVVtYQNAEru3/Z/py5j+H/3Yt3eHGkCJac1a9YsdOvWDc8995zUoSiarCptuSmd0dC1pM7HGOrYaMq2P23CJFXCYcnzho7i4O26sIMnERGRfuv25mD+ttMoLq8y+lhBAIrLKzF/m+ZylVE9Wtk5OlKCAwcOYO3atTh+3HjXdbKO7CptRERERGSdo7mFmL/tjEkJm67i8irM33YGx/IK7RQZKUVZWRnGjh2LJUuWoEWLFlKHo3iyStqCh5xARYUrKipqz2fL3tgF2Ru7GDzWlG1/UnajdGaWbOsE6k9XxYzVMchYHSN1GERERKIvdmeipKLSomNLKiqxfHemjSMipVm4cCFCQkIwYsQIqUOpF2S1PRLQlOe1dAdttx5x1KbPo71WjNsPjbM00a0vXRXDRh+SOgQiIiLRtaJS/JZ+tdp7KnMIAvDr2au4XlTKrpKk16lTp/DZZ5/h8OHDUKlUUodTL8iq0kZERERE1kk5lGf1OVQAUg5bfx5SnsrKSowdOxb//Oc/ERwcLHU49YbsKm1lJW7i162GmVfhOZ8cAQBoOdT4xZCWVtgyVsfUWVnJ29wJQYNP1j7ORnPaiIiIiOpy5tLtWl0izVVSUYUffj+CZleOoEmTJmjSpAm8vb3Fr728vFhhUbC6ZvqtW/EV3N3dMWHCBKnDrFdkl7QRERERkeVul1TY5DwXr93ExsPbcevWrVr/KisrxQRO959uYmfon/YxjRo1YuInM8Zm+i35+SzuZRbiy48+h4sLN+w5kuySNrW7ZQtNdlIkWg+17Norc65vM3b9kr4qG2D+9V0Xt3QEADz8/CmzjiMiIqL6zdvTNm/vHu/xCJYuHa/3vtLSUjGBu337tt7E7tKlSzh79qzB+0tLS9G4cWOjyV1d/xo3bszkwUY0IyLOoKSiUu/1kCX3Ezh16xi8s+Mibrl6czSEA8kuadNtepG1XpNMmdoZ0hjdxia65NKMRHeOnG6ylrn2flKpM1TblLl0REREVP+09/eGh/qSVVskPdUuaB/Q2OD9Hh4e8PPzg5+fn8XPUV5eLiZ0hhK7GzduIDs72+Bj7t69i4YNG5pU2avrfrVadm+JHcqcmX6AijP9JFC//0KJiIiIFGZITBCW7ky36hwCgCFdg2wTkAFubm5o1qwZmjVrZvE5KisrcefOnToTv1u3biE/P9/gY+7cuQNPT0+ztnbq++fu7m7DV8dxrJ3pFxnkg8ggHztFR1r1Kmkzt7GJI5xPjhAbpxiqFupW2LRYZSMiIiJ9mjfywBPhLfDz6csWtf1XqYDe7Vo4Rbt/V1dX+Pj4wMfH8qRBEAQUFRUZ3e6ZkZFR5/1qtdrs6/pq/vP09HT4dX62mOn31ahYG0dFNck6adOXmGRtiELbkWl2f+7spEiHDOI2pdOlHLEbpvJd2KTpxhrygnP+jRIR1Wev9grF7xnXUFxu/ptxT7UrJvcKtUNU8qRSqdC4cWM0btwYQUGWVRcFQUBxcbHRxO/8+fN1Pqaqqsqipi66/xo2bGhy4seZfs5D1kkbEREREZmvS7APZj7T3ozrlDS83Fww85n23O5mJpVKhQYNGqBBgwYICAiw+Dy6DV4MJXYFBQU4c+aMwftLS0trJXOGEsC00haoqvS07meHZqbfK4+3teo8VDerkralS5fi66+/hkqlQkREBFatWoWCggIMHz4c169fR0xMDNauXWvTPb5ClQrXvg8HADR/zrr92nVxRJXNGRiqODprlS3/204IHKS/wyegaVajbxutoyqvcuLsFTYp1iciIjnRNoioqyOglkqlqbDNfKY9G0tIyNYNXgwldtevX8e5c+dw0L0zyhu1sSrmkooqnCm4Y9U5yDiVIFhWEM3Pz0fPnj1x6tQpeHl5YejQoXjmmWewbds2DBo0CMOHD8fEiRPRpUsXTJo0qc5zeat80V3Vp9ptGWu6IuylB4mBoc6Pzs5R2z2JHG2fsAu3hRuSPLe91ycicm5Srk9SOJZXiOW7M/Hr2atQ4UHrdkDTJVKA5hq2yb1CWWGrZ8YkHsAvZ65YfZ4+7f2wIv4RG0REhlhVaauoqEBxcTHc3Nxw7949BAQE4JdffsGGDRsAAPHx8Zg9e7bRN0VERLbG9YmISCMyyAdfjYrF9aJSpBzOw5mCO7hdUg5vTze0D2iMIV2DeD1SPWWrmX7enm42OQ8ZZvFvKjAwEG+99RZCQkLg5eWFvn37IiYmBj4+PuKsi6CgIOTn55t13rzNnQAAYYOrb79zVIVNd1aaI0hdZdO+3oaGghM5I3utT0REzqxZIw9ed0TVOGKmH9mGxSPkb968ia1btyI7OxsXL17E3bt3sX37dpOPT0hIQGxsLGJjY1GOUvH24iIPFBdJ92lP6+HH6s21S9lJkQgafJIJGymOvdYnIiIiJRkSY/0sPkfM9CMrkradO3eidevWaNGiBdzc3DBo0CD8+eefKCwsREVFBQAgLy8PgYGBeo+fMGECDh48iIMHD8INLMkTke1wfSIiIjJOO9PP0tFwzjTTz9lZnLSFhIRg7969uHfvHgRBwK5du9CxY0f07t0bKSkpAIDExEQMGDDArPO6eVbAzbOi2m2Za6MtDdNseZs7iVsGla6+VBSp/rHX+kRERKQ0r/YKhafa1aJj69tMPylZnLR1794dQ4YMQdeuXREREYGqqipMmDABCxcuxJIlSxAaGorr169j7NixZp3X1bUKrq5VyNoQZWloZsv/thPyv31wbZczbxfMWu+4BJdIruy1PhERESmNdqafl5t5aQFn+jmWxS3/bUm3pXZuSmcAQFmZWmzSkbk2GqEvHrHb82sTtrrmdzmLrPXRaBtnv9eKyFRKaanNlv9EyqOU9YnIltbtzTFpph8gwMtNzZl+Dia7pE2frPXRcHPXbJmsrHBB6xFHARgflExE0lHKmyImbUTKo5T1icjWjM30q6yqQsm5Q0iZOx7dwwKkC7Qess1wBjvJXKfZ6ldV7CpWj7I3dhHvNyVhO3d/m2UbDrAmIiIiIjLIlJl+k8Yk4rcta9H9nXekDrdekXXSRkREREREjlXXTL9Zs2ahd+/emDRpEho35nw2R5F10hY66sG1WWLFbIR5FTPdCpu2uUldA621jTzkcl1YdlKkTbs8XtgUgZAXjtvsfERERERUf3Tq1AlPPvkkPv/8c7z33ntSh1NvWNw90l4Mtdx39yiHu0e5VeduOzKtzoQNAFxcK+HiWmnxc+hu36wpY01XZKzpatb5bN2WnwkbEREREVlj1qxZWLp0KW7duiV1KPWG7JI2IiIiIiKSr3bt2uGZZ57BsmXLpA6l3pBd90ht8xGhSoWwlw5LGZbJzidHoOXQ2hWsvM2dLJr5lr2xi9ghk+Qta0OU3uptzjeRAIBWw5Q9wFw7+F7fSA6ldGdj90gi5VHK+kQkpczMTPTo0QMZGRlo2rSp1OEoHittRERERERkltDQUAwYMABLliyROpR6QXaVNmNyUzojeMgJvfel//cRAIBX02KDj7G3nG8iFV9dkZqhyibJi1I+yWaljUh5lLI+EUktOzsbsbGxSE9PR7NmzaQOR9GcLmmzVNaGKAhVKgDVu1LKyZWt7QEAfgPOmH1s+spYAED4mIM2jYmcg627jNqCUt4UMWkjUh6lrE9EcjBx4kT4+PhgwYIFUoeiaNweSUREREREFnnvvfeQkJCAK1euSB2Kojldpa2u5h51NUWwhLZ9P5uCSI9bIp2PUj7JZqWNSHmUsj4RycVrr70GT09PLF68WOpQFMspkrb8bzshcJD5XRi1xwKw+Hhnd2FTBADOZ7OEHLccOhOlvCli0kakPEpZn4jkIj8/HxERETh58iQCAgKkDkeRuD2SiIiIiIgsFhgYiPj4eF7XZkdOkbRZUyULHHRSPD43pTNyUzrbKiyjtNsrpSQImn+2cG5DlG1O5CSUXGU7nxwhdQhERESkINOnT8fatWuRl5cndSiKJLvtkdrEoI3OwGJnGjZdkNoBAQNPm/z4ukYYEDkzpWw/4vZIIuVRyvpEJDdvv/027t27hy+++ELqUBTHKSptREREREQkb++88w6SkpJw4cIFqUNRHFkkbaWtvZC1XtP5sc3ItGpVNqB698bspEiHxmYuc6psAFhlIyIiIiJFaNGiBV555RXMmzdP6lAURxZJm0d2MdrGmdamX6rrjLI3drHZNWrGrqvLWB1jk+ep+Zz6njd9Zaw4mJvIGuc2RNW76x6JiIioujfffBObN2/GuXPnpA5FUWSRtBERERERkfNr1qwZXn31VVbbbEx2jUi01SBuGyS5KUjtAMD8LbD1lVIu9GcjEiLlUcr6RCRXN2/eRFhYGPbs2YOwsDCpw1EEtdQB1GQoWdPd2seEjqTAZI2IiIjIuKZNm2Lq1KmYO3cu1q5dK3U4imB0e+SYMWPg5+eHzp0fJE03btzAU089hbCwMDz11FO4efMmAEAQBLz++usIDQ1FZGQkDh8+bL/Iiaje4/pEREQkT1OnTsX27dtx+jQ/9LYFo0nb6NGjsX379mq3LViwAH369EFGRgb69OkjTj//6aefkJGRgYyMDCQkJGDSpElmB5S5LhqZ66Jr3R485IT4zxQXNkXgwibTBghnrI6xWfMPNmIgchxHr09ERERkGm9vb/zjH//A3LlzpQ5FEYwmbY8//jh8fX2r3bZ161bEoM1qEgAAH8ZJREFUx8cDAOLj45Gamire/tJLL0GlUqFHjx4oLCxEQUGBeREJmn/WtvYPeeE4Ql44blJCFjb6EMJGH7Lq+bRqjisgIvtx+PpEREREJpsyZQp++eUXnDjBS5usZVH3yMuXLyMgIAAA4O/vj8uXLwMA8vPzERwcLD4uKCgI+fn5es+RkJCA2NhYxMbGohylloRBRFQL1yciIiJ5aNSoEd5++23Mnj1b6lCcntUt/1UqFVQqldnHTZgwAQcPHsTBgwfhBg/x9tAXjyD0xSMG57Fd3NLRrOcxVkWz1ew1Miz/205Sh0D1lK3XJyIiIjLP5MmT8eeffyItjbvRrGFR0vbQQw+J24oKCgrg5+cHAAgMDERubq74uLy8PAQGBtogzAcefv6UTc/XesTRat9nro1G5tra19TVJcvM69jOJ5t2rZ096Lte0N4CB510+HNS/SXl+kRERETVNWjQANOnT2e1zUoWJW39+/dHYmIiACAxMREDBgwQb1+zZg0EQcDevXvRpEkTcZsSEZEjcH0iIiKSl1deeQUHDhzAoUO26SFRHxkdrj1ixAjs3r0b165dw0MPPYQ5c+Zg4MCBGDp0KC5cuICWLVsiOTkZvr6+EAQBr732GrZv344GDRpg1apViI2NNRqEkofXartJGmpQcim1A/zvz//SVsFCRx1xTHBEduSI4bVcn4jIEhyuTeR4X3zxBbZt24Yff/xR6lCcktGkzRGU/KaISRvVV0p5U6Tk9YmovlLK+kTkTEpLSxEWFobk5GT06NFD6nCcjtMkbdq2/Sp1FZMaM2ivn2s59LjEkVB9o5Q3RUzaiJRHKesTkbP5z3/+g2+//RY7duyQOhSnI+ukzdaVJ22nSC+vMrG6BQA532hmwrUapr9jJREAsUFN6Iv80MAUSnlTxKSNSHmUsj4ROZuysjK0a9cOa9euRc+ePaUOx6lY3fKfiIiIiIjIGHd3d3zwwQf48MMPpQ7F6ciu0pa3WTPTK2iw4TbxuSmdAQDBQxw3Xd3YtWm6tCMA2prw2Pog55tIVjHrIaV8ks1KG5HyKGV9InJG5eXl6NChA/773/+id+/eUofjNNRSB1BTaalbrduyNkShqswVgGZYtrFkzR5JginJmhaTteqYsBERERERALi5ueHDDz/ErFmz0KtXL6hUKqlDcgrcHklERERERA4zcuRIXL58Gbt27ZI6FKchq6TtwqYItB2ZVqtS1XZkGsJGH0LY6ENiM4i6OHtlR9vx0d6ykyId8jz2kJ0UKZv4c76JFJvZZG2IErfHEhEREVFtarUas2fPxgcffAAZXKnlFGR3TVv6fx4BAIS/ckD2XR2luLZOH14zRgCQ/rVmUHT4uIMSR6KhlGtGeE0bkfIoZX0icmaVlZWIjIzE4sWL8fTTT0sdjuzJqtJGRERERETK5+rqitmzZ+PDDz9ktc0EskvaVF6VUHlVImNNV7QadkzWFaTgISckrbJpt+XJ+TUCnHsbpjMJH3dQNlU2IiIiImMGDx6MsrIy/PDDD1KHInuy2x5pzIVNEQh54bidIyJnV5DaAQE6A9SN4eBs21PK9iNujyRSHqWsT0RKkJqaijlz5uDQoUNwcZFdPUk2+MoQEREREZEkBgwYAJVKhdTUVKlDkTXZVdrSv+wGAAiftN/g47PWa6oibeNYFSGSK6V8ks1KG5HyKGV9IlKKH374ATNmzMDRo0dZbTNAfq+KWtD8qyFjdQwyVscA0CRr9TFh0yarJB/narT31/07JSIiIiLj+vXrhwYNGiAlJUXqUGRLfkkbERERERHVGyqVCnPnzsXs2bNRWVkpdTiypJY6gJrCxx/Qe7tHwzIHRyI/9bK6eL+SVXPgutS0FbY2NeIKG31IinCIiIiInFrfvn3RtGlTJCUlIS4uTupwZMdpKm0hLxxn10gJZCR21X+7g7YAth2ZJlnCph3urk+bkWm1EjYiIiIisoy22jZnzhxUVFRIHY7sOE3SRkREREREyvW3v/0NDz/8MNavXy91KLIju6Qtc220ODOLzJdVozGGtcLiD+u/vR5sA5T70HIiIiIiJVGpVJgzZw7mzp2L8vJyqcORFVkkbWVtvHA+OQKAZrixdsDxpdQOuJTaQcrQnI7crv1yZnUlwOwSSURERGR7TzzxBFq3bo3ExESpQ5EVWSRtREREREREADB37lzMmzcPZWVsRKhlNGkbM2YM/Pz80LlzZ/G2t99+G+3bt0dkZCSef/55FBYWivd9/PHHCA0NRbt27bBjxw6TgnA/V4yWQ48jOymyWgXDf+Bp+A88be7PRAqS801knQ1B7KmuqmXY6EP1Youo3DlifSIiIiLHeuyxx9ChQwesXLlS6lBkw2jSNnr0aGzfvr3abU899RROnDiBY8eOITw8HB9//DEA4NSpU0hKSsLJkyexfft2TJ48mbMWiMhuuD4REREp05w5czB//nyUlJRIHYosGE3aHn/8cfj6+la7rW/fvlCrNSPeevTogby8PADA1q1bMXz4cHh4eKB169YIDQ3F/v37TQ6m9fBjeisYmeseNCbJ3tjF5PPVlLe5E/I2d7L4eHK8VsOOsSEIGeTI9YmIiIgcp1u3boiKisJ///tfqUORBauvaVu5ciWefvppAEB+fj6Cg4PF+4KCgpCfn2/tUyB01IOh0mq3SmStj0bWevM7TAYNPomgwSetjseZZCdFIjtJmu2FRFJzxPpERERE9jFnzhx8/PHHKC4uljoUyVmVtM2fPx9qtdqiqeUJCQmIjY1FbGwsylFqTRhERLVwfSIiInJuXbt2RY8ePfDVV19JHYrk1JYeuHr1avzwww/YtWsXVCoVACAwMBC5ubniY/Ly8hAYGKj3+AkTJmDChAkAAG+Vr97H6FN80wvh4w9YGrbJ0lfFIPxl52800Xo4txZS/SPV+kRERES2NXv2bPTt2xcTJkxAw4YNpQ5HMhZV2rZv345Fixbhu+++Q4MGDcTb+/fvj6SkJJSWliI7OxsZGRno1q2bzYIF4JCEDYCkCduOi5y1pnRPnrgjfr05by++yz+A7/Id87etdFKuT/UB16f6hesTEUktMjISjz/+OL744gupQ5GU0UrbiBEjsHv3bly7dg1BQUHi3tLS0lI89dRTACCWLTt16oShQ4eiY8eOUKvV+OKLL+Dq6mr3H4KI6ieuT0RERMo3a9Ys9O7dG5MmTULjxo2lDkcSKkEQBKmD8Fb5oruqj9RhENmdtsK2s3NjsRNq6xFHcT45AgDQcuhxi84759whzGoTY5sgbWSfsAu3hRtSh2E1rk9EyqOU9YmoPomLi0OnTp3w3nvvSR2KJJi0EUkge2MXtB5xFAAQdQRIM78Zql2dT46wOIHUUsqbIq5PRMqjlPWJqD45e/YsevbsiczMTDRp0kTqcBzO6pb/RERERERE9tSuXTs888wzWLZsmdShSMLi7pFEZL7NeXsBAIODNBU2QH5VNqD2Nk3tXMS2cUf0PZyIFODB+tRD4kiIiPT74IMP0KNHD7z++uto2rSp1OE4FCttRA7kpnKFm8oV55MjkBYtz4RNn7ZxR5iwESmcdn0iIpKr0NBQDBw4EEuWLJE6FIdj0kZERERERE7h/fffx/Lly3Ht2jWpQ3EoJm1EDtQ/8BH0D3zE5CYf6atikL7qQVdIbcfJmrTbF4mILKVdn4iI5KxVq1Z44YUXsHjxYqlDcSgmbUQyVZDaAWqPSqg9KnEptQMupXaAm3uF3seaunUxOykS2UmRZseSvbFLrYTxwqYIs89DREREZK2ZM2ciISEBV65ckToUh2HSRkRERERETiM4OBhxcXFYtGiR1KE4DLtH0v/f3t0HNXXmewD/BtEo11IHlZcmVAiR1xBAWPQvXbbV7dpZ6ktrQVzr0tau3k472tG161bA0XWvne7Y2Tq9w3bW+gp2vVthKr7ccW1rd8pFwJcuviRiUnlTQWWobhWQ5/6RJhVNCAnJOYf4/fwjOTmcfDnB3/DkOef3kMK07U8CAETNOYfGPen9ntPOb3D6PRd3ZmB0SDcAoKc7GH19ts9jHpyBs28H0G9x7we1fpqMJ+aedTx2ts+TLwxtHTciIiIib7399tswGAx46623EBUVJXccv+Pi2kQKYtqWiWD1PQCAABC38NSQjtfy9xRo5jkf6LljHzDGLTzluGcueGTvoO/HC5TFa1mfiAJPoNQnokfdihUr0NfXh/fff1/uKH7HyyOJiIiIiGjY+e1vf4udO3eiublZ7ih+x5k2IoW58sPlkQAQOedcv+ese42IefHMQ99jKTciNu/h7QBwcZdtlky/yHmzkvsvx/z2E1tzkcHOpg0kUD7JZn0iCjyBUp+ICFi9ejVu376NrVu3yh3FrzjTRkREREREw9KqVatQXl6Oy5cvyx3FrzjTRqQglrI0R1v/+5uOWPfa2vTfP8tmKUtz2iBEKo170ge85y5QPslmfSIKPIFSn4jI5ne/+x06OjpQWloqdxS/4aCNiBzaKxMAABNzL7jdl4M2IhquAqU+EZHN9evXER8fjxMnTkCn08kdxy94eSQREREREQ1b48ePx+uvv44NGzbIHcVvOGgjGkbsa6s9tL3c6JPjT8y9MKhZNmDoyxEQERER+cqKFStQWVkJs9ksdxS/4KCNaBhxdQ/b6NE9/R5b9xod98HZmbZlDnhse+fIB/lqQEhERETkL+PGjcObb76J9evXyx3FLzhoIyIiIiKiYe/NN9/E4cOHce7cOfc7DzPBcgcgGiz7jI+r9ciGk5JLdSjSPTzz1bg7A3EFtvXULu7MgP5Xtq+d/eyWcqNjhi1qzrl+HSbV6h9n3uzNRXBN9ePr7EmH6LM9dvzb1f8zHPt6cbFzzvTbZu9cFHXfGnKX9qRDx8sliYiISEahoaFYuXIl1q9fj7KyMrnj+BS7RxINwPRRFgAg/pVayV+7+X9S+rX9H4irRbflFCjd2VifiAJPoNQnInrYrVu3EBcXh6NHj8JgMMgdx2d4eSQREREREQWEsWPHYtWqVSguLpY7ik+5HbQVFhYiPDzc6Uj1vffeg0qlQkdHBwBACIE33ngDer0eRqMR9fX1vk9MJKH4V2plmWUDgJ7uwV+9LMcsmxIalLA+ERER0YOWL1+Of/7znzh1KnBu3XA7aFuyZAkOHTr00PampiYcOXIETz75pGPbwYMHYTabYTabUVpaimXLlvk2LdEjpK/P/US4vUtk2w/3n0lJCfcWsj4RERHRg0JCQrBmzZqAmm1z+1fh9OnTERYW9tD2FStWYPPmzVCpfmxuUFFRgcWLF0OlUmHatGno7OxEW1ubbxMTEf2A9YmIiIicWbp0KWpra1FXVyd3FJ/w6p62iooKaDQapKX1X+i3paUF0dHRjsdarRYtLS1Oj1FaWoqsrCxkZWWhB3e9iUEU0OIKTqLl7ylo+XuKy31iXjyDmBfP9OvkKBXzxwOv+yYX1iciIiIaM2YM3n77baxbt07uKD7h8aDt3//+N/7whz8MeeE6++i3trYWI6Ee0rHsmvYZ0LQvcLrEEGnmNUAzrwEXd2Xg4q4Ml/s52vpLaPIS5X1ypeT6RERERNJ65ZVX8M0336C6ulruKEPm8aCtsbERFosFaWlpiImJQXNzM6ZMmYIrV65Ao9GgqanJsW9zczM0Go1PAxMRucL6RERERHZqtRq///3vUVRUJHeUIfN40Jaamopr167BarXCarVCq9Wivr4ekZGRyM3NxY4dOyCEQHV1NR5//HFERUX5I7dT0c//C9HP/0uy1yOSin7RSegXnXxou2lbJkzbMnHz2mMypFIeJdcnIiIikt6SJUtgMpnw1VdfyR1lSNz2FM/Pz8fnn3+Ojo4OaLValJSU4OWXX3a67+zZs1FVVQW9Xo+QkBBs27bN54GJAp29lX5fXxDiFg7cqjb+19Jcomi/NNPZwFFOrE9EREQ0kFGjRuGdd97BunXr8I9//EPuOF5TCSGE3CFCVWGYqnpK7hhEiuDJoE0q3gza/k8cRZe44a9IkmF9Igo8gVKfiGhwenp6kJSUhL/85S/IycmRO45XvOoeSaQU336SKncEn4vNO4PYvDMeDdga96T7MZHryzOJiIiIlG7kyJEoKipCUVERFDBf5RUO2mhYm7TgG7kj+I2lLM39Tj8QfSr3OxERERE9ovLz83Ht2jUcPXpU7ihe4aCNiIiIiIgCWnBwMIqKivDOO+8My9k2DtqIFCo2/7Tj67b9SWjbn+RyX860EREREQ1swYIF6OrqwqFDh+SO4jEO2oiIiIiIKOCNGDECJSUlWLdu3bCbbeOgjWgYiJpzDlFzzj20/dtPUvHtJ6kQ3fyvTEREROTOvHnz0N3djc8++0zuKB5xu04bEcnP3iVz0oJv0F6ZAACYmHtBskYs1r22ZQhiXjwjyesRERER+UNQUJBjtu3ZZ59FUNDw+OB7eKQkIiIiIiLygeeeew5BQUHYv3+/3FEGjYtrE8nk/tkzZ1o/TcYTc89KGcljjbtti27HFTy8hlugLF7L+kQUeAKlPhGR9w4cOIA1a9bg9OnTw2K2TfkJadi6/LdUXP5b4C1+7SuTFnwz4OWNT8w9i8Y96YNaOPvKAJ0l/Smu4KTTARsRERGRks2ePRshISHYt2+f3FEGhYM2IiIiIiJ6pKhUKqxfvx7FxcW4d++e3HHcYiMS8psnX5CmSUYgi1t4asDnLeW2BiGxc/zbIMS0LRMAEP/rOr++DhEREZFUZs2ahbCwMJSXl6OgoEDuOAPiTBuRzOz3hdlZytJgKUtzPPfg8/eLzTuD2LzBD9gu7srAxV2uj2fda3R0irQzbctE/K/rOGAjIiKigGKfbSspKUFvb6/ccQbEmTYiIiIiInok5eTk4IknnsDu3bvx7Pw87KtrxvkrXei604vQ0cFIjAzFC5lajB+rljWnortHmj7MBgDEL6uROhKRIrjrMGl3ZX8SIp0svi2nQOnOxu6RRIEnUOoTEfnGx5XHUPy3rzEqJgMqAHd7+xzPjQ4OggDw04SJWD5Dj7TocbJkVPRMGwdr9Kgb7OLZsn/yQkRERDQM7aq24r9q7wKaVHTfN1izu/PDtiNnr+JLUwfWzk7EomkxEqdU+KCNiIiIiIjIH3ZVW7Gx6hy+7+kDVAO3+hAC+L7nHjZW2a5sknrgxkYkCtQm05pbpAwDrW9nX7ftwbXbou67NNLeURIAzB9n+ickERER0TB2uqkTG6vO2wZsHvi+pw8bq87jTHOnn5I5p+h72ogu7rR1OtT/igs4S6lxd8aQF80OlHtGWJ+IAk+g1Cci8t7SnbX433NX4c1ISKUCfp4cgf9elOX7YC5wpo2IiIiIiB4ZHbfu4gtTu1cDNsB2qeSxC+24fuuub4MNgIM2UjT9r05ylm0QLj1wueRQDXWWjYiIiEip9tU1D/kYKgD76od+nMFyO2grLCxEeHg4DAZDv+1//vOfkZiYiJSUFKxevdqxfdOmTdDr9UhISMDhw4d9n5hoAJf2pPt8ADMc6BaekjuCLFifiIiIyFPnr3T1a+vvjTu9fTjf9p2PErnntnvkkiVL8Prrr2Px4sWObceOHUNFRQVOnz4NtVqNa9euAQDOnj2L8vJyNDQ0oLW1FU8//TRMJhNGjBjhv5+AiB5ZrE9ERETkqa47vT46To9PjjMYbmfapk+fjrCwsH7bPvzwQ6xZswZqtW1l8PDwcABARUUF8vLyoFarERsbC71ej5oa5ay1Zt1rhHWv0f2Oj6hAOD+6hace2VmnR1Eg1SciIiKSRuho36x6Fjp6pE+OMxhe3dNmMplw/PhxTJ06FTNmzMCJEycAAC0tLYiOjnbsp9Vq0dLS4pukPhDz4hnEvHhG7hgw75gidwSnlHJ+iIZiuNYnIiIikkZiZCjUwUNr7TE6OAiJUY/5KJF7XqXt7e3FjRs3UF1djXfffRcLFiyApysHlJaWIisrC1lZWeiBdJ1XiCiwsT4RERHRQJ7P1A75GALA81OGfpzB8mpuUKvVYt68eVCpVMjOzkZQUBA6Ojqg0WjQ1NTk2K+5uRkajcbpMZYuXYqlS5cCACZMmADTf9Rg4sSJ3sTxq/b2dp/nmrDFgptDnGzzRy5fUWo25vLMUHONssrTnNbX9Wns2LG4mWiRJLunAvV3x1+Yy3NKzTZc6xMRKcOEsWrMiJ84pHXachImYvxYte/DueDVoG3OnDk4duwYcnJyYDKZ0N3djQkTJiA3NxcLFy7EypUr0draCrPZjOzsbLfH6+joQFZWFmpra72J41fM5TmlZmMuzyg1lzu+rk+JiYmKPQ9KfY+YyzNKzQUoN5tScxHR8PGfP9XjuLkD3/fc8/h7RwePwPKf6v2QyjW3g7b8/Hx8/vnn6OjogFarRUlJCQoLC1FYWAiDwYBRo0Zh+/btUKlUSElJwYIFC5CcnIzg4GBs3bqVndmIyG9Yn4iIiMgbadHjsHZ2IjZWncP3PYNv/z9mZBDWzk6EUTvOj+ke5nbQVlZW5nT7rl27nG5fu3Yt1q5dO7RURESDwPpERERE3lo0LQYAsLHqPO703hvwUkmVyjbDtnZ2ouP7pDSiuLi4WPJXdSEzM1PuCE4xl+eUmo25PKPUXFJT8nlQajbm8oxScwHKzabUXEQ0vBi14zB98gTcvN2NppvfY2SQCr19P47eRgcHYUSQCk8nhWPzfCNmJkfKklMlPG2rRkREREREFGCu37qLffXNON/2Hbru9CB09EgkRj2G56doJW064gwHbURERERERAome8/bQ4cOISEhAXq9Hn/84x9lzdLU1IScnBwkJycjJSUF77//PgCguLgYGo0G6enpSE9PR1VVleTZYmJikJqaivT0dGRlZQEAbty4gZkzZ2Ly5MmYOXMmbt68KWmmCxcuOM5Jeno6QkNDsWXLFtnOV2FhIcLDw2EwGBzbXJ0jIQTeeOMN6PV6GI1G1NfXS5pr1apVSExMhNFoxNy5c9HZ2QkAsFqtGDNmjOPc/eY3v5E010Dv3aZNm6DX65GQkIDDhw/7LZeSKKU+sTZ5hrXJ+1ysTURECiVk1NvbK3Q6nWhsbBR3794VRqNRNDQ0yJantbVV1NXVCSGE6OrqEpMnTxYNDQ2iqKhIvPvuu7LlEkKISZMmifb29n7bVq1aJTZt2iSEEGLTpk1i9erVckQTQtjey4iICGG1WmU7X1988YWoq6sTKSkpjm2uztGBAwfEM888I/r6+sTXX38tsrOzJc11+PBh0dPTI4QQYvXq1Y5cFoul337+5CyXq/euoaFBGI1GcefOHXHp0iWh0+lEb2+vJDnloqT6xNrkPdYmz3KxNhERKZOsM201NTXQ6/XQ6XQYNWoU8vLyUFFRIVueqKgoTJliW/X6scceQ1JSElpaWmTL405FRQVeeuklAMBLL72E/fv3y5bl6NGjiIuLw6RJk2TLMH36dISFhfXb5uocVVRUYPHixVCpVJg2bRo6OzvR1tYmWa5Zs2YhONjWvHXatGlobm72y2t7msuViooK5OXlQa1WIzY2Fnq9HjU1NX5OKC8l1SfWJu+xNnmWi7WJiEiZZB20tbS0IDo62vFYq9Uq5g8Rq9WKkydPYurUqQCADz74AEajEYWFhZJf6gMAKpUKs2bNQmZmJkpLSwEAV69eRVRUFAAgMjISV69elTyXXXl5OfLz8x2P5T5fdq7OkZJ+9/7617/iF7/4heOxxWJBRkYGZsyYgePHj0uex9l7p6TzJRWl/sysTZ5hbfIeaxMRkXLIfk+bEt26dQvz58/Hli1bEBoaimXLlqGxsRGnTp1CVFQU3nrrLckzffXVV6ivr8fBgwexdetWfPnll/2eV6lUUKlUkucCgO7ublRWVuKFF14AAEWcL2fkPEeubNy4EcHBwSgoKABgm1G5fPkyTp48iT/96U9YuHAhurq6JMuj1PeObFibPMPa5D3WJiIiZZF10KbRaNDU1OR43NzcDI1GI2MioKenB/Pnz0dBQQHmzZsHAIiIiMCIESMQFBSEV199VZZLL+znJTw8HHPnzkVNTQ0iIiIcl820tbUhPDxc8lwAcPDgQUyZMgUREREAlHG+7FydIyX87n388cf47LPPsHv3bscfbGq1GuPHjwdgW4MoLi4OJpNJskyu3jslnC+pKe1nZm3yHGuTd1ibiIiUR9ZB209+8hOYzWZYLBZ0d3ejvLwcubm5suURQuDll19GUlISVq5c6dh+//0En376ab+OVlK4ffs2vvvuO8fXR44cgcFgQG5uLrZv3w4A2L59O5577jlJc9mVlZX1u/xI7vN1P1fnKDc3Fzt27IAQAtXV1Xj88ccdlypJ4dChQ9i8eTMqKysREhLi2N7e3o579+4BAC5dugSz2QydTidZLlfvXW5uLsrLy3H37l1YLBaYzWZkZ2dLlksOSqpPrE3eYW3yHGsTEZFCydsHxdYpa/LkyUKn04kNGzbImuX48eMCgEhNTRVpaWkiLS1NHDhwQCxatEgYDAaRmpoqfvnLX4rW1lZJczU2Ngqj0SiMRqNITk52nKeOjg7xs5/9TOj1evHUU0+J69evS5pLCCFu3bolwsLCRGdnp2ObXOcrLy9PREZGiuDgYKHRaMRHH33k8hz19fWJ5cuXC51OJwwGgzhx4oSkueLi4oRWq3X8nr322mtCCCH27dsnkpOTRVpamsjIyBCVlZWS5hrovduwYYPQ6XQiPj5eVFVV+S2XkiilPrE2eY61ybtcrE1ERMrExbWJiIiIiIgUjI1IiIiIiIiIFIyDNiIiIiIiIgXjoI2IiIiIiEjBOGgjIiIiIiJSMA7aiIiIiIiIFIyDNiIiIiIiIgXjoI2IiIiIiEjBOGgjIiIiIiJSsP8HPoduT00v+fgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1080x360 with 3 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "for m in masks: \n",
    "\n",
    "    fig, (ax1, ax2, ax3) = plt.subplots(1, 3, figsize=(15,5))\n",
    "    plt.title(str(m))\n",
    "    \n",
    "    # Full adjacency\n",
    "    ax1.set_title('Full Adjacency mask')\n",
    "    adj = show_adjacency_full(m, ax=ax1)\n",
    "    \n",
    "    # Filtered adjacency\n",
    "    f = filter_adj(adj)\n",
    "    ax2.set_title('Full Adjacency mask')\n",
    "    ax2.imshow(f);\n",
    "    \n",
    "    # Plot subgraph\n",
    "    ax3.set_title(\"Subgraph\")\n",
    "    G = nx.from_numpy_array(f)\n",
    "    G.remove_nodes_from(list(nx.isolates(G)))\n",
    "    nx.draw(G, ax=ax3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
